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Abstract: Resurgent transseries have recently been shown to be a very powerful construction in 
order to completely describe nonperturbative phenomena in both matrix models and topological 
or minimal strings. These solutions encode the full nonperturbative content of a given gauge or 
string theory, where resurgence relates every (generalized) multi-instanton sector to each other 
via large-order analysis. The Stokes phase is the adequate gauge theory phase where an 't Hooft 
large expansion exists and where resurgent transseries are most simply constructed. This paper 
addresses the nonperturbative study of Stokes phases associated to multi-cuts solutions of generic 
matrix models, constructing nonperturbative solutions for their free energies and exploring the 
asymptotic large-order behavior around distinct multi-instanton sectors. Explicit formulae are 
presented for the Z2 symmetric two-cuts set-up, addressing the cases of the quartic matrix model 
in its two-cuts Stokes phase; the "triple" Penner potential which yields four-point correlation 
functions in the AGT framework; and the Painleve II equation describing minimal superstrings. 
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1. Introduction and Summary 

For almost 40 years the large N limit of gauge theory [1] has been source to many fascinating 
results and ideas, with large duality [2] playing a definite central role. Given some nonabelian 
gauge theoretic system, this limit produces an asymptotic perturbative expansion^, in powers of 
I/N"^ (or, from the point -of-view of large duality, in powers of the closed string coupling g^). 
What this means is that the genus g perturbative contributions to the free energy^, Fg(t), will 
display large-order behaviour of the type Fg ~ (2(7)! and the perturbative expansion will have zero 
radius of convergence [3]. But, more importantly, what this implies is that the perturbative series 
is not enough to define the free energy and nonperturbative corrections of the type ~ exp {—N) 
are needed in order to properly make sense out of this expansion. 

The study of nonperturbative corrections, their relation to the large-order growth of the 
perturbative expansion and their application within the resummation of perturbation theory has 
a long history. It is also almost 40 years since these topics were first considered within the 
quantum mechanical context of the quartic anharmonic oscillator [4] . Later, they were extended 
to the usual perturbative expansion in quantum field theoretic systems; see, e.g., [5] for a review 
of early developments. In the present work we are interested in yet another extension of these 
ideas, namely towards trying to understand, nonperturbatively, the expansion. This topic 
has a more recent history where we will follow previous work in [6, 7, 8, 9, 10, 11, 12, 13, 14, 15] 
(we refer the reader to the introduction of [15] for a quick overview of these developments, 
or to [16] for an excellent review of some of the main ideas which were put forward in the 
aforementioned references). For the moment let us simply mention that most of the original 
ideas and results in [5] have their counterparts within the large context. In particular, the 
leading growth of the free energies Fg(t) is dictated by a suitable instanton action, A{t), whose 
physical origin is associated to eigenvalue tunnelling, at least within the matrix model context 
[17, 18, 7]. Moreover, the subleading growth is associated to the one-loop amplitude around the 
one-instanton sector [17, 18, 7]. Further corrections to the large-order {2g)\ growth arise from 
higher loop amplitudes around some fixed instanton sector — corrections in 1/ g — and from higher 
instanton numbers — corrections in l/n^, where n is the instanton number [15]. 

A very interesting novelty emerges as one addresses the asymptotic perturbative expansion 
around some fixed multi-instanton sector. Although one might naively expect that the large- 
order growth of this asymptotic expansion would be controlled by sectors with different instanton 
numbers (either higher or lower), it turns out that this expectation is incomplete: one actually 
needs to introduce new nonperturbative sectors in order to match all large-order results [13, 15]. 
In the examples studied so far, these "generalized" multi-instantons have instanton actions with 
opposite sign as compared to the "physical" multi-instantons, and they may be assembled all 
together into a so-called transseries solution to the problem at hand — where resurgent analysis 
relates every generalized nonperturbative sector to each other via large-order analysis. In this 
way, by checking that no further corrections exist deep in the large-order asymptotics, one is led 
to the assumption that these resurgent transseries solutions completely encode the full nonper- 
turbative content of a given gauge or string theory. Transseries solutions were first introduced in 
the string theoretic context in [8], and this approach was later extended in [15] (see, e.g., [19] for 
a mathematical overview). The use of resurgent analysis as a tool towards fully understanding 
generalized nonperturbative sectors of string theoretic systems was first introduced in [20, 13], 

^This is, of course, the very well known topological or genus expansion due to 't Hooft [1]. 
^If one has the gauge theory in mind, in here t = Ngs is the 't Hooft coupling. If one has the dual closed string 
theory in mind, in here f is a geometric modulus associated to the background geometry. 
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and those approaches were later extended in [15] (see, e.g., [21, 22] for mathematical overviews). 
In particular, we refer the reader to [15] for a complete exposition of the ideas and techniques we 
shall use in the course of the present work. In fact, it is precisely one of our main goals in this 
paper to extend the techniques of resurgent analysis and transseries solutions to other examples 
beyond the ones in [15]. Furthermore, let us mention that ideas from resurgence and transseries 
have also appeared in, e.g., [23, 24, 25], in a quantum mechanical context, and have recently 
been shown to be very promising in the study of quantum field theory [26, 27, 28, 29, 30, 31]. 

In our endeavor to understand the nonperturbative structure of the large limit, we have 
began with a simpler class of gauge theoretic systems: random matrix models^. In this case, it 
is very well known that, at large N , matrix eigenvalues cluster into the cuts of a corresponding 
spectral geometry [32]. Within the context of large duality, this was later understood as 
relating matrix models to B-model topological string theory in local Calabi-Yau geometries 
[33, 34, 35] (see [36] for a excellent review). Depending on the potential appearing in the matrix 
partition function, and the phase in which the model is to be found, generically the spectral 
geometry will correspond to a multi-cut configuration. This is an important point, specially in 
light of the question: is it always the case that, as one considers the large N limit of some given 
gauge theory, one will find an expansion of 't Hooft type with a closed string dual? Within the 
matrix model context, this question was first raised in [37] and answered negatively. 

Let us dwell on this point for a moment as it is also at the basis of the class of examples we 
choose to address in this work. The nature of the large N asymptotic limit depends very much 
on which gauge theoretic phase one considers [37, 12] and is analogous to the study of Stokes 
phenomena in classical analysis. When considering single cut models, one finds the familiar l/N"^ 
expansion [38, 39, 40], i.e., one finds a topological genus expansion with a closed string dual. 
However, this is not usually the case when considering multi-cut models, where one finds large 
theta-function asymptotics instead [37, 41, 10], i.e., there is no genus expansion and possibly 
no closed string dual. These two distinct large N asymptotics are associated to what we call 
Stokes and anti-Stokes phases, respectively, generalizing the usual Stokes and anti-Stokes lines 
in classical analysis [16]. In fact, in the Stokes phase the single-cut would correspond to the 
leading saddle, with pinched-cuts corresponding to exponentially suppressed saddles [7]. On 
the other hand, in the anti-Stokes phase the many cuts correspond to many different saddles 
of similar order, where their joint contribution translates into an oscillatory large N behavior 
[37]. Of course one should start the analysis in the opposite direction: having identified Stokes 
and anti-Stokes phases with particular large N asymptotics, one may then ask what spectral 
geometry configurations appear in each distinct phase. The point of interest to us in here is 
that there are regions of moduli space where the Stokes phase is actually realized by a multi-cut 
configuration (essentially, by configurations where all cuts are equal, i.e., they have the precise 
same eigenvalue filling). It is this type of multi-cut configurations which we investigate and 
explore in this work, within the framework of resurgent transseries. 

This paper is organized as follows. We begin in section 2 by reviewing background material 
concerning both matrix models and resurgent transseries. We briefly review saddle-point and 
orthogonal polynomial approaches to solving (multi-cut) matrix models, and then go through the 
basics of resurgence and transseries. In section 3 we address the multi-instanton analysis when 
we have two cuts with Z2 symmetry (ensuring we are in a multi-cut Stokes phase). This is done 
using methods of spectral geometry which essentially generalize previous work in [7, 9]. Elliptic 
functions and theta functions which, due to the elliptic nature of the spectral curve, appear 

■^Although further motivated by their relation to topological strings, these models are probably the simplest 
already showing all features of large A'^ resurgence, without the added complication of renormalon physics [27, 28] . 



-3- 



during the calculation, end up canceling in the final result thus providing further evidence on 
the nature of the Stokes phase. This is actually an interesting point of the calculation, as, on 
what concerns the perturbative sector, it was source to some confusion in early studies of Z2 
symmetric spectral configurations. In fact, the original saddle-point calculation of the two-point 
resolvent in a Z2 symmetric distribution of eigenvalues [42, 43], with an explicit elliptic function 
dependence, did not match the corresponding orthogonal polynomial calculation [44, 45, 46], 
which saw no trace of these elliptic functions. The reason for this was that [42, 43] worked in a 
fixed canonical ensemble, while in the spectral curve approach to solving some given multi-cuts 
scenario one needs to address the full grand-canonical ensemble as later explained in [37]. We 
shall explicitly see in section 3 what is the counterpart of those ideas within the multi-instanton 
context. Section 4 presents applications of these general multi-instanton results into two different 
examples. On the one hand we study the two-cuts phase of the quartic matrix model. This 
example further explores the quartic matrix model along the lines of [15], in particular as we 
construct a two-parameters transseries solution in this phase. Instantons in this example are 
associated to S-cycle eigenvalue tunneling [7] and we provide tests of our analytical results by 
comparing against the large-order behavior of perturbation theory. Earlier results addressing the 
asymptotics of this model were presented in [47] and we extend them in here within the context 
of transseries and resurgent analysis. On the other hand, we address the example of the "triple" 
Penner potential which is associated to the computation of four-point correlation functions within 
the AGT framework [48]. This is an interesting example as it is actually exactly solvable via 
generalized Gegenbauer polynomials and where instantons are associated to 74-cycle eigenvalue 
tunneling [11]. In section 5 we turn to the asymptotics of multi-instanton sectors, and we do this 
in the natural double-scaling limit of the two-cuts quartic matrix model, which is the Painleve II 
equation. This equation describes 2d supergravity, or type OB string theory [49, 50, 51, 52], and 
we fully construct its two-parameters transseries solution, checking the existence of generalized 
multi-instanton sectors via resurgent analysis. Earlier results addressing the asymptotics of this 
model were presented in [53] and we extend them in here within the context of transseries and 
resurgent analysis. In particular, we compute many new Stokes constants for this system (in this 
way verifying and generalizing the one known Stokes constant [54]), and present the complete 
nonperturbative free energy of type OB string theory. We close in section 6 with a discussion of 
some ideas which could lead to future research. Let us also stress that, due to the nature of the 
large-order analysis, we have generated a large amount of data concerning both the two-cuts 
quartic matrix model and the Painleve II equation. Due to space constraints it is impossible to 
list all such results in the paper, but we do present some of this data in a few appendices. 



2. Revisiting Multi-Cuts Matrix Models 

Let us begin by setting our notation concerning both saddle-point and orthogonal polynomial 
approaches to solving matrix models, with emphasis on multi-cut configurations. We shall also 
review the required background in order to address the construction of (large N) resurgent 
transseries solutions for these multi-cut configurations, when in their Stokes phases. 

2.1 The Saddle Point Analysis 

Let us first address the saddle-point approximation to computing the one-matrix model partition 
function (within the hermitian ensemble, /3 = 1) in a general multi-cuts set-up; see, e.g., [32, 55, 
42, 36, 9]. In such configurations the eigenvalues condense into s different cuts Ci U • • • UCg = 
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[xi,X2] U • • • U [x2s-i, X2s], and, in diagonal gauge, the partition function is written as 



Z{Ni,...,Ns) 



Ni\---Nsl 



A. •••/() n(|^)^^(^^)^"-^-^'^''^ (2-1) 



with 't Hooft couphng t = Ngs (fixed in the 't Hooft Umit). In the above expression the {A^^^} 
are the eigenvalues sitting on the I-th cut, with kj = 1, . . . , Nf and Yl'i=i — A(Aj) is 

the Vandermonde determinant. In this picture it is natural to consider the hyperelliptic Riemann 
surface which corresponds to a double-sheet covering of the complex plane, C, with precisely the 
above cuts. One can then define ^-cycles as the cycles around each cut, whereas i?~cycles go 
from the endpoint of each cut to infinity on one of the two sheets and back again on the other. 
For shortness, we shall refer to C as the contour encircling all the cuts, i.e., C = IJ/=i 

The large N saddle-point solution is usefully encoded in the planar resolvent, defined in 
closed form as 

' 2tJc 27ri z-wV asiw)' ^ ' 



where we have defined 



2s 



(Ts{z) = \{{z-Xk) (2.3) 

k=l 

and where one still needs to specify the endpoints of the s cuts, {xk}- One may now describe 
the large A'^ matrix model geometry via the corresponding spectral curve, y{z), which is given in 
terms of the resolvent by 

y{z) = V'{z)-2tujo{z) = M{z)^/^Jz). (2.4) 

If the potential V{z) in the matrix model partition function (2.1) is such that V'{z) is a rational 
function with simple poles at z = /3j, i = 1,2,..., A; and with residues ai at each pole, the 
expression for M{z) in the expression above is simply 

r d^n^ J— +y — (2.5) 

At this stage one still needs to specify the endpoints of the cuts. If the eigenvalue distribution 
across all cuts is properly normalized, the planar resolvent will have the asymptotic behavior 
^o{z) ~ ^ as z — 7- +00. In turn, this asymptotic condition implies the following set of constraints 

-"^^^^ = 21 Sr.., (2.6) 
27ri y/as{w) 

with n = 0,l,...,s. These are s + 1 conditions for 2s unknowns, where the remaining s — 1 
conditions still need to be specified and they come from the number of eigenvalues Nj one chooses 
to place at each cut. This distribution of eigenvalues may be equivalently described by the partial 
't Hooft moduli tj = QsNj, which may be written directly in terms of the spectral curve: 

ti = ^^j^Azy{z), I=l,2,...,s. (2.7) 
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Notice that, as expected, these are only s — 1 conditions as they are not ah independent, i.e., 
"^j^i tj = t. Both constraints (2.6) and moduli (2.7) now determine the fuU spectral geometry. 
It is also useful to define the holomorphic effective potential 

VLd^) = yiz). (2.8) 

In this case, the effective potential is given by the real part of the holomorphic effective potential, 
in such a way that 

yeff(A) =Me / dzyiz). (2.9) 



2.2 The Approach via Orthogonal Polynomials 

While saddle-point analysis is the appropriate framework to describe the spectral geometry 
of multi-cut configurations, it gets a bit more cumbersome when one wishes to address the 
computation of the full free energy. In the 't Hooft limit, where N — )• +oo with t = QsN held 
fixed, the perturbative, large A^, topological expansion of the free energy is given by^ 

+00 

F{gs,{ti}) = logZ^ Y.9l'~^^9iti)- (2-10) 

9=0 

Computing this genus expansion out of a hyperelliptic spectral curve has a long history — starting 
in [38, 42], passing through [37], and recently culminating in the recursive procedure of [39] — and 
it is in fact an intricate problem in algebraic geometry [40]. 

An easier approach to computing the free energy of a matrix model is to use the method of 
orthogonal polynomials; see, e.g., [56, 55, 36, 15]. On the other hand, this method is less general 
as it is not applicable to arbitrary multi-cut configurations. However, as we shall also see in the 
course of this paper, orthogonal polynomials do work when addressing multi-cuts Stokes phases. 
As such, let us swiftly review this method in the context of the one-cut solution (the multi-cuts 
extension will be addressed later). Considering the partition function (2.1) with a single cut, one 
may consider the positive-definite measure on M given by 

d^(z) =e-^^(^)^. (2.11) 

Normalized orthogonal polynomials with respect to this measure are introduced as Pn{z) = 
+ ■ ■ ■ , with inner product 

dfj.{z) Pn{z)pm{z) = KSnm, n > 0. (2.12) 

As the Vandermonde determinant may be written A(Aj) = detpj_i(Aj), the partition function 
of our matrix model may be computed as 

N~l N 

Z=l[K = hl^l[r^--, (2.13) 

n=0 n=l 

where we have defined r„ = , ^" for n > 1. These r„ coefficients further appear in the recursion 
relations 

Pn+l{z) = {z + Sn)Pniz) - rnPn-l{z), (2.14) 



'Throughout this paper we shall use the symbol ~ to signal when in the presence of an asymptotic series [15]. 
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together with coefficients {sn} which will vanish for an even potential. Plugging the above (2.14) 
in the inner product (2.12) one obtains a recursion relation directly for the r„ coefficients [56]. 

One example of great interest to us in the present work is that of the quartic potential 
V{z) = + ^z^. In this case it follows that Sn = and [56] 



A 
6 



-1 + Tn + r„+i 



ngs- 



(2.15) 



The free energy of the quartic matrix model (normalized against the Gaussian weight Vq{z) 

1 , 

2' 



hz^, as usual) then follows straight from the definition of the partition function (2.13) 



+ 0O 



N 



^ . F - ^ log I . Y./r^T,it) 4 log I . I i E (1 

^ g=0 ^ U ^'^ n=l 



n 
N 



log 



(2.16) 



This genus expansion is made explicit by first understanding the large N expansion of the r„ 
recursion coefficients. Changing variables as x = ngs, with x G [0,t] in the 't Hooft limit, and 
defining TZ{x) = r„ with TZ^{x) = x, the above example of the quartic potential (2.15) becomes 
[56, 15] 

(2.17) 



n{x) + ^ {T^{x - Qs) + n{x) + n{x + g,)) | 



x. 



As TZ{x) is even in the string coupling, it admits the usual asymptotic large A'^ expansion 

+ 00 

n{x)c^Y.sl'R2g{x), 
9=0 



(2.18) 



allowing for a recursive solution for the R2g{x). In particular, in the continuum limit the sum 
in (2.16) may be computed via the Euler-MacLaurin formula (with i?2fe the Bernoulli numbers 
and X = t£,) 



N 



n=l 



1 



yielding 



Hi. 9s) 



t 

2is 



2N 



ho TZ{x) 



5=1 +00 



1 B2k 



-0 ^A^^M2fc)! 



^(2fc-l)(^) 



+00 

9=1 



23-2 -P2ff 



(2^)! dx^9~i 



+ 

x=o/ 9. 
{t — x) loj 



1 /■* 

-7T / dxit — x) log 
li Jo 



n{x) 



?=i 



5=0 



+ 



(2.19) 



n{x) 



x=t 



x=0 



(2.20) 



This analysis was first presented in [56] and was recently extended to a full resurgent transseries 
analysis in [15], and we refer the reader to these references for further details. We shall later see 
how it generalizes to accommodate the two-cuts Stokes phase of the quartic matrix model. 



2.3 Transseries and Resurgence: Basic Formulae 

The discussion so far has focused upon the large A, perturbative construction of the matrix model 
free energy (2.10). If, on the other hand, one wishes to go beyond the perturbative analysis in 
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order to build a fully nonperturbative solution to a given matrix model, one needs to make use 
of resurgent transseries. This subject was recently thoroughly addressed in [15], and we refer the 
reader to this reference for full details on these techniques and their origins. In here, we shall 
nonetheless cover just enough background to make the present paper a bit more self-contained. 

Resurgent transseries essentially encode the full (generalized) multi-instanton content of a 
given non-linear system and, as such, yield nonperturbative solutions to these problems as expan- 
sions in both powers of the coupling constant and the (generalized) multi-instanton number (s). 
In general, many distinct instanton actions may appear and, as such, transseries will depend 
upon as many free parameters^ as there are distinct instanton actions. For most of this paper, 
and similarly to what was found in [15] for the one-cut quartic matrix model and the Painleve 
I equation [13], a two-parameters transseries will be sufficient to describe the two-cuts quartic 
matrix model and the Painleve II equation. These two-parameters transseries generalize the 
one-parameter cases which were first introduced in the matrix model context in [8]. 

Similarly to what happened in [15], we shall only need to consider the special case of a 
two-parameters transseries ansatz with instanton action A and "generalized instanton" action 
—A. This may be written as 

+ 00 +00 

F{z,ai,a2) = E E <^i<F^''^"'Hz), (2-21) 

n=0 m=0 

where z is the coupling parameter (here chosen ~ l/gs) and cri, a2 are the transseries parameters. 
Further, the above (n|m) sectors label generalized multi-instanton contributions of the form 

+00 p(n|m) 

F(-l-)(^) ^ ^-in-m)Az $(^|^^^(^) ^ ^-in-m)Az ^ (2.22) 

9=1 



In this expression f3nm. is a characteristic exponent, to which we shall later return when needed. 
Resurgent transseries are defined along wedges in the complex z-plane (upon Borel resummation, 
see, e.g., [15] for details) and they are "glued" along Stokes lines in order to construct the full 
analytic solution. This "gluing" is achieved via the Stokes automorphism &g which essentially 
acts upon the transseries (2.21) by shifting its parameters. For instance, given an one-parameter 
transseries with Stokes line on the positive real axis, the gluing is achieved by shifting a — )■ a + Si 
as one crosses from the upper to the lower positive-half-plane, where Si is the Stokes constant 
associated to that particular Stokes line — although, generically, there may be an infinite number 
of distinct Stokes constants. In our two-parameters case, there are two sets of Stokes coefficients, 

(k) ~(k) 

and , labelled by integers k and i with A; > 0. Do notice that not all of these coefficients 
are independent and in [15] some empirical relations between them have been found, in the 
Painleve I context. We refer the reader to that reference for further details. 

The main point of interest to us in this subsection concerns large-order analysis [5] , and how 
resurgent analysis improves it [13, 15]. Recall that if a given function F{z) has a branch-cut in 
the complex plane along some direction 9, being analytic elsewhere, then 

1 r'"'"" Disc gF(w) 
Fiz) = / dw ^ , (2.23 



''Free parameters which are essentially parameterizing the corresponding nonperturbative ambiguities. 
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where we have assumed that there is no contribution arising from infinity [5, 4]. The key point 
now is that the aforementioned Stokes automorphism &g, which may be expressed as a muhi~ 
instanton expansion [15], relates to this branch-cut discontinuity as 

6e = l-Disc0, (2.24) 

in such a way that the discontinuity itself may be written in terms of multi-instanton data. For 
instance, starting with the perturbative sector, it was shown in [15] that in the two-parameters 
transseries set-up (2.21) there will be two branch-cuts, along 6 = and 6 = n, such that with 
/3oo = one finds 

+ 00 



Disco$(o|o) = -Y.(^n e-'^'^W0)> (2-25) 
Disc^$(o|o) = -E (2-26) 



k=l 

+ 00 

E 

k=l 

Using (2.22) and (2.23), we then find the perturbative asymptotic coefficients to be given by [15] 

k 



77(0|0) _ V ' J ^ - ^k,0) ^ ^ \.y - MM -n^ij (fc|o) ^^ ^ 

(^-i) r (g - /3o,fc) r (g - /3o,fc -h + i) (o|fc) 

What this expression shows is that the asymptotic coefficients of the perturbative sector, for large 
g, are precisely controlled by the coefficients of the (generalized) multi-instanton sectors, (n|0) 
and (0|n). Of course that besides the coefficients F^"'^^ and Fg^^^\ the perturbative coefficients 
also depend on the two Stokes constants, S^^ and S^l, and these still need to be determined. 
For the moment, let us just note that the leading large-order growth is dictated by the Stokes 
constants and the one-loop (generalized) one-instanton coefficients F^^^^^ and F^^^^\ Higher loop 
coefficients in the (1|0) and (0|1) sectors yield corrections in 1/g, whereas the higher (n|0) and 
(0|n) sectors yield corrections which are suppressed as 

As we turn to the models of interest to us in the present work — such as matrix models 
or topological strings — there are a few extra points to consider. First, the perturbative sector 
(2.10) is given by a topological genus expansion, in powers of 2g — 2, where the string coupling 
is z = l/gs- Secondly, as we address matrix models or topological strings, one needs to consider 
a version of the multi-instanton sectors (2.22) where both the action A and the perturbative 
coefficients Fg"'^'^^ become functions of the partial 't Hooft moduli (or geometric moduli) tj. 
But, more importantly, due to resonance effects which will later appear in either the quartic 
matrix model or the Painleve II equation, one also needs to consider the inclusion of logarithmic 
sectors as [13, 15]: 

A(tr) ^""^ r,, Adr) 

k=0 g=0 

(2.28) 
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nm 



We shall later uncover that the maximum logarithmic power is k'j^m — ^mn — min^Ti, ttx) — in (5^ 
and that /3nm = /3mL = /3(m + n) — [{knm + ^)/2]/, where [•]/ is the integer part of the argument 
and where /3 = 1/2. In practice, this essentially means that all we have done up to now was for 
the k = "sector" , and that the /3nm coefficients take into account the fact that the perturbative 
expansions may in fact begin at some negative integer. Going back to the perturbative (0|0) 
sector in (2.10), we know that F^"!^) is given by a genus expansion containing only powers of the 
closed^ string coupling g'^. Thus, one needs to impose i^2^+i = in (2.27), which will produce a 
series of relations between the (0|A;) and (A;|0) contributions since its right-hand-side must vanish 
order by order in both ^ and k~^. As further explained in [15], in the end we find that for all k 
and g, 

(5f ) = (-l)^+<i (2.29) 

When working out the full details of either the two-cuts quartic matrix model or the Painleve 
II equation, we shall find further relations between different coefficients iT^"'"*^!*^] ^ either when m 
and n are exchanged, or relating the k ^ coefficients to the k = coefficients. In some cases, 
these will imply further relations between different Stokes constants. 

Finally, using the above relations (2.29) back in the large-order formula for the perturbative 
sector (2.27), we obtain the asymptotic large-order behavior of the perturbative coefficients in 
the string genus expansion (2.10) as 



This procedure may be extended in order to find the large-order behavior of all (generalized) 
multi-instanton sectors. In particular, we are here interested in the large-order behavior of the 
physical multi-instanton series 

i?("|o), xhe precise calculation is a bit more cumbersome due to 
the logarithmic sectors appearing in (2.28), and we refer the reader to [15] for full details. The 
final result is 

f('^|o)[o] _ + ^^ (Sj^^)'' r(g + /gj^^o - /gj+fc^o) ^ r(g + - Q - h) ^n+k\om (,..h 

^ ~ ^ V n / 27ri /, ^ t-/ , ^ ^ ^ ^ 

fc^i V n / zTTi (^kAY+f^nfi Pn+kfi T{g + p'^'^Q - /3;,|;t,o) 

+00 I ^ k 



k=i y m=i e=o 7ier(m,fc) 5jer{m,m-e+i) \j=i 

kn+t-k,l -pf . n[0] o[r] \ +00 




E''^* '^'^ ^(g + /^n,0 Pn+l-k/) r(g + /3„^Q Pn+l-k/ ^) j^(n+i-k\£)[r] ( i A\h 

r=o i-kAY+^-:° Pn+e-k,i r(5 + /3|^;o - /3^|^_;,^^) 



6ro + Q{r-l)(BkA{a) + daf ^ BkA{a) 



(2.31) 



Let us define the many ingredients in this expression (but, again, we refer the reader to [15] for 
the full details). The sums over 7j and 5j are sums over Young diagrams, where a diagram 7^ G 



One has to be a bit careful with the precise meaning of the labels: in full rigor, the coefficients Fg in (2.10) 

,(0|l 

29 



precisely stand for -^29'"' in the present transseries language, as can be seen by comparing against (2.22) 
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T{k,i) : < 7i < • • • < 7fc = ^ has length i{T) = k, and where the maximum number of boxes 
for each ji is ^(r^) = £. The sum over 6s is similar, now with 0<6i<62<---<6k = k — m + 1 
and < 5s < s + 1. These 6s form a diagram T(k,k — m + 1) that has length £(T) = k and 
^(r-'") = k — m + 1, with an extra condition that each component 6s G r(A;, k — m + 1) has at 
most s + 1 boxes. For these definitions to be consistent we still have to set 70 = 0, (^o = 1) a-s well 
as the Stokes constants S^^ = S^^ = S^^l = 0. Next, defining d7j = 7^ —7^-1, and similarly for 
d(5s, one has 

S(n,i) = ((j + 1 - 6,) Si'^j;) + (n - 7,. + J + 1 - 6,) S^^^^'^^) © (j + 1 " '^.) , (2-32) 

where G(x) is the Heaviside step-function. Finally, we have introduced the function 

Bs{a) = ^(a + 1) - log(-s) = Bs{a) - (2.33) 

with il^i^) = digamma function. 

There are a few relevant features to be found in (2.31). Besides the multitude of (generalized) 
multi-instanton sectors and Stokes constants that now play a role, there is also a new type of 
large-order effect. In fact, and unlike the usual perturbative case which had a leading large-order 
growth of gl, essentially arising from the Gamma function dependence, we now find a large-order 
growth of the type g\ log arising from the digamma function, and this is actually a leading effect 
as compared to the g\ growth. The first signs of this effect were found in [13], in the context of 
the Painleve I equation, and further studied in [15]. 

3. Multi— Instanton Analysis for Z2— Symmetric Systems 

Having reviewed the main background ingredients required for our analysis, we may now proceed 
with our main goal and address the nonperturbative study of Stokes phases associated to multi- 
cuts configurations. These phases arise when all cuts are equally filled and, to be very concrete 
and present fully explicit formulae, we shall next focus on two-cuts set-ups (see [9] as well). In 
this case, equal filling also implies that the configuration is Z2-symmetric. As we shall see in 
detail throughout this section, this symmetry implies that hyperelliptic integrals which appear in 
the calculation will reduce to elliptic integrals, and that, physically, the system will be found in a 
Stokes phase. Notice that, strictly within the orthogonal polynomial framework, it was already 
noticed in [57] that equal filling of the cuts would lead to a Stokes phase. 

3.1 Computing the Multi Instanton Sectors 

Let us begin by considering the multi-instanton sectors of a two-cuts matrix model. We shall do 
this by following the strategy in [9], i.e., we shall consider the two-cuts spectral geometry as a 
degeneration from a three-cuts configuration. In principle one could also consider degenerations 
from more complicated configurations if one were to introduce several distinct instanton actions, 
but for our purposes degenerations from three-cuts will suffice. In this case, a reference filling of 
eigenvalues across the cuts is of the form ( A^i , N2 , N3 ) , with A'^i + A^2 + -^3 = the degeneration 
will simply be N2 — ?• 0, and the Z2 symmetry will eventually demand A^i = N3. 

In matrix models, (multi) instantons are associated to (multiple) eigenvalue tunneling [17, 
18, 7, 9] and, as such, the multi-instanton sectors are described by tunneling eigenvalues in- 
between the three cuts, as shown in figure 1 (it is simple to see that two integers, ni and 77-2, are 
enough to parameterize all possible exchanges of eigenvalues between three cuts, i.e., all possible 
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Figure 1: Eigenvalue tunneling as the multi-instanton sectors of a three-cuts matrix model. 



background choices). In the particular case of the Z2-symmetric two-cuts configuration, the 
reference background is of the form 

'N N\ , , 

-,0,- . (3.1) 



As we shall see later on, the one-instanton sector will correspond to summing over all configura- 
tions which leave a single eigenvalue on the middle-cut. From the spectral geometry viewpoint, 
the Z2 symmetry essentially places the cuts at [—6, — a]U[a, b] and the spectral curve (2.4) becomes 

y{z) = M{z)y^{z^ -a^){z^ -b"^), (3.2) 

where M{z) is given by (2.5). In this configuration, the pinching cycle will be found at z = 0. The 
action associated to eigenvalue tunneling essentially measures their energy difference in-between 
cuts [17, 18, 7, 9], as given by the holomorphic effective potential (2.8), and in the particular case 
of this Z2-symmetric configuration with equal filling it is simple to check that the equal filling 
essentially translates to 

/■a 

dzy{z) = 0. (3.3) 

This condition will further imply that one may completely evaluate all data in the spectral 
geometry just by using the asymptotics of the resolvent (2.6). One is left with one instanton 
action to evaluate, describing tunneling from each of the (equal) cuts up to the pinched cycle^ 
located at xq such that M{xq) = [7]. In here xq = and 



A= rdzy{z). (3.4) 

J a 

Having briefly explained the set-up, one may proceed and compute the partition functions 
associated to the relevant configurations along the lines in [9]. Let now y{z) be the spectral 
curve (2.4) of the three-cnts configuration, with cuts [a;i,X2] U [x3,X4] U [a;5,a;6]- Let us consider 
the aforementioned set-up with Ni — ni, N2 + {ni + 712) and N-^ — 722 eigenvalues in the first. 



^This is the non-trivial saddle located outside the cut, where eigenvalues may tunnel to [7], 



- 12 - 



second and third cuts, respectively, and let us consider the associated multi-instanton amplitude 
written in terms of the 't Hooft moduli (2.7) 



^■(ni.na) = ^ {h - mg^, ^2 + Uig^ + ra2ffs, ^3 - n2gs) 

Z (tl, ^2, ^3) 

with ti + t2 + = t. For convenience we introduce the variables 

si = ^ (ti -t2- h) , (3.6) 

S2 = ^ (ts -t2- h) , (3.7) 

and use them to expand the exponent of (3.5) above (i.e., the difference of free energies between 
the "eigenvalue-shifted" configuration and the reference background), around gs = and for 
^1,112 <^ N. One simply finds* 



= exp (^-^Y^mdsTo^ exp riin^dsA^Foj |l + 0(ff.)}. 



(3.8) 



In this expression we find two, in general different, actions 

Ai = ds,Fo, i = l,2, (3.9) 

which may be computed in terms of geometric data if we use the special geometry relations 

dFo 



dzy{z). (3.10) 



In the present three-cuts configuration, the two actions are then given by 



dFn 

M = ^= / dzy{z), (3.11) 

1x2 



dsi 

(■Xi 

ds2 



A2 = ^= I dzy{z), (3.12) 
' X5 



and they have the usual geometric interpretation appearing in figure 2, generalizing the one-cut 
case appearing in [7, 9]. The extension to an arbitrary number of cuts is straightforward. The 
other feature we find in (3.8) are the second derivatives of Fq, and for those it is convenient to 
introduce the (symmetric) period matrix 

= ^ 3.13 
27ri osiOSj 

Having understood the general form of the multi-instanton amplitudes, we still need to 
understand the precise nature of the multi-instanton expansion. The grand-canonical partition 
function is obtained as a sum over all possible eigenvalue distributions into the multiple cuts, 



'^For shortness we shall many times omit the arguments; it should be clear that whenever we write Fq we always 
mean the reference configuration -Fb(ii, ^2, fs), and similarly in other cases. 
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Figure 2: The two-cuts spectral curve y{z) is a genus one curve, with a pinched cycle at the non-trivial 
saddle xq which is obtained by taking X3 — > 0:4 in the three-cuts geometry (in the Z2 symmetric scenario 
Xq — 0). The instanton actions Ai and A2 naturally appear as i?-cycles in this spectral geometry. 



with their total number fixed. In our case, and making use of the multi-instanton amplitudes 
(3.5), this translates to 

Z{N) = Yj Y1 ^^"^'"'^ (3.14) 

nt=—N2+N3 n2=—ni—N2 

Let us now consider the reference background of interest to us, i.e., the Z2-symmetric two-cuts 
configuration describing a multi-cuts Stokes phase. This background has moduli ti = = ^ 
and t2 = 0, in which case both instanton action will be equal Ai = A2 = A, as well as th = T22. 
Changing variables from ni and n2 to i = ni + 71-2 and m = ni — 712, we may write the multi- 
instanton amplitudes (3.8) as 

= exp (^-y^ exp (^-^ (rn + r^) exp (m - r^) m^^ |l + 0(g,)|, (3.15) 

where it now becomes clear that it is ^ = ni + ^2 > which will label the multi-instanton 
sectors. Of course this further implies that we still need to sum over the "relative" index m in 
order to obtain the "purely" ^-instanton amplitude: it is this sum over m which essentially moves 
our calculation to the grand-canonical ensemble. In other words, the grand-canonical partition 
function (3.14) is of the schematic form 

ZiN) = ^(^=0) + Z(^=i) + ^(^=2) + . . . = Z(^=°) (^1 + IS + IS + • • • ) ' (3-^6) 

where now each term Z^^^ contains a sum over all possible values of m = ni — n2 that satisfy 
ni + n2 = i. Fixing i eigenvalues on the middle-cut implies that we only have N — i available 
eigenvalues to place in each of the two side-cuts, which yields the limits on the m-sum. But 
because m jumps by values of two, it turns out that it is actually more convenient to change 
variables and use as the "relative" index m = 2r — i. Overall, we find 



(3.17) 
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With a certain abuse of notation, we shall immediately identify the i~th. instanton amplitude as 



where all that is now missing is the explicit evaluation of the many different ingredients which 
appear above, in particular explicitly evaluating the sum. 

Let us begin by addressing the period matrix (3.13), i.e., the second derivatives of the planar 
free energy. Using the special geometry relation (3.10) and the explicit form of the spectral curve 
(2.4), it follows that 

|^ = (-iy«r"d.(-2)M, (3.19) 

dsidsj J^^^ dsi 

where the derivative of the resolvent has the form^ 

dituoiz)) _ Ci;\t,sk) + C['\t,sk)z 



dsi yJ(T-i{z) 



(3.20) 



The coefficients which appear in this expression, CQ\t, s^) and c[^\t, s^), may be fixed by taking 
derivatives of the partial 't Hooft moduli (2.7), and by using the definition of the variables {si}, 
(3.6) and (3.7), as 




dz -" .^^^2, 1=1,2,3. (3.21) 

9si \_^J 2-k\Jai V^3(^) 

Note that although this relation corresponds to a system of 6 equations for 4 unknowns, two of 
the equations are redundant as we can deform contours in order to find = (there is no 

residue at infinity) . If we now define the integrals 

K,i= I — — ^ and Ci= I> — — (3.22) 
J A' 2vri ^a-i{z) J a' 27ri ^^a2,{z) 

(i) 

then we can express all the coefficients Cj in terms of these integrals as 

(1) _ C1+C2 (2) _ ^2+^3 no\ 

LiK,2 — i-21^1 1-31^2 — i-21^3 

(1) _ /Ci +/C2 (2) _ ^2 +/C3 . . 

L'21^1 — l-\f^2 L-2f^3 — L'31^2 

So far these results are only formal: hyperelliptic integrals are hard to evaluate. However, they 
may in fact be explicitly evaluated when one imposes Z2 symmetry into the problem. In this 
case, one places the cuts as [—6, —a] U [— c, c] U [a, h] (where we shall later be interested in the 
c — )• degeneration) and it immediately follows that 

/Ci = /C3 = -^/C2 = -/C, (3.25) 

^In order to check this relation one exphcitly uses (2.4) and (2.5) when taking derivatives, and this will yield 
the polynomial structure in z. In order to fix the degree of this polynomial, one compares the asymptotics as 
z — >■ +00 on both sides of the equation. Generically, the degree will depend on the number of cuts as s — 2. 
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= -£3 = -C, £2 = 0, (3.26) 

leading to the (simplified) coefficients 

^0'^ = Ci>^ = ^' (3-27) 
-C['^ = ) = ^. (3.28) 
As they will be needed in the following, let us also introduce the B-cycle integrals: 

r.r^=r^, (3.29) 

J-a ^/<73{z) Jc ^/<T3{z) 

C^- rdz-^= Pdz-^ (3.30) 

All these A and S-cycle integrals may be explicitly evaluated, and expressed in terms of complete 
elliptic integrals of the first kind, K{k'^), with k being the elliptic modulus. This is also the 
technical reason why one may find Stokes phases within multi-cuts configurations: symmetries 
(in this case a Z2 symmetry) may effectively reduce hyperelliptic geometries to elliptic ones! The 
results are 



dx 



,2 



^ J-b ^ VWsixYi ttVo:^^ ^ w - ay ' ^^-^^^ 

and, for (3.29) and (3.30), 

1 fb^c^-a^)\ ^ ^ 

C = , ^ K ( • (3-34) 

Having explicitly evaluated all integrals, we may now start assembling these results back 
into our original formulae and address the degeneration limit c — )• 0. In order to do that, it is 
first important to notice that this limit must be taken carefully as the free energy is not analytic 
in the 't Hooft modulus associated to the shrinking cycle [9]. This may be explicitly seen by 
splitting the free energies as 

Fg{h,t2,h) = F^{t2) + Fg{ti,t2^3), (3.35) 

where F^[t2) are the genus g free energies of the Gaussian model depending on the vanishing 
't Hooft modulus, which, at genus 5 = and g = 1, have a dependence as logi2- As explained 
in [9], for the ^-instanton sector it is not appropriate to look at the integration over the i 
eigenvalues in the collapsing cycle as a large approximation; rather one should exactly evaluate 
the Gaussian partition function associated to this cycle, which is 
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with G2 {i + 1) the Barnes function. Then, the partition function around the £-instanton con- 
figuration should be properly written as 

zW = ZfZ^'\ (3.37) 

where all "hatted" quantities in Z^^"^ are now regularized and analytic in the t2 ^ limit. 

The instanton action is the simplest quantity to evaluate as it is in fact regular in the t2 ^ 
limit. One simply finds 

A= r dzM{z)^{z'^ - a'^){z^ - b'^){z^ - c'^) > r dz M{z)y^{z^ - a^){z^ - b^), (3.38) 

Ja '^^O Ja 

where A4{z) = zM[z). To compute the period matrix we must first address the second derivatives 
of the planar free energy, (3.19), which are given by 

,2 ^ _ .2 ^ 1 r dz 1 r , z ic c 



and by 



aiFo . a;.Fo = ^ + ^ ^ = - . -, (3.39, 



8.,8„Fo = |-^. (3.40) 

With these results, the period matrix follows immediately. In particular we obtain 

i 

rii+ri2 = -, (3.41) 

i 

rii-ri2 = ---. (3.42) 

TT L 

The need for regulation of the shrinking cycle is now very clean. In fact, if one takes the c — )■ 
limit in (3.41) above one obtains 

lim ivr (m + T12) ~ 2 lim log c + log H . (3.43) 

c^O c->0 \ 160^62 J 



However, as explained, this logarithmic divergence — which emerges in one of the elliptic integrals — 
will be precisely cancelled by the "Gaussian divergence" arising from the shrinking cycle. The 
regulation is simply [9] 

d^Jo = lim (afFo - log t2) , (3.44) 
where the vanishing 't Hooft modulus is, via (2.7), 

1 



t2 = — dzM{z)^/{z^ - a2)(z2 _ 62)(^2 _ c2). (3.45) 

Changing variables z — ?• z/c, expanding the result in powers of c and performing the integration, 
it follows^o 

t2 = ^Mabc'^ {l + 0{c^)) , (3.46) 
'^'^In the purely three-cuts scenario it is simple to check that M is just a constant; more on this in the following. 
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which win indeed cancel the divergence above. As for the combination (3.41), it has a regular 
c — 7- limit. Using known properties of elliptic integrals [58] one may compute 

i K{1 -e) _ iK' . 

where the elliptic modulus in this Z2-symmetric limit is simply given by A; = 

Finally, in order to obtain the multi-instanton amplitudes (3.18), all one has to do is evaluate 
the sums in (3.17). When £ = 0, the sum in (3.17) yields the Jacobi (elliptic) theta-function 
given by 

+00 

d2,{z\q) = 1 + 2^ q"'^ cos (2rz) . (3.48) 

r=l 

In fact, using this definition it is straightforward to evaluate 



N/2 

^hm^ E exp(|^(Tn-Ti2)(2r)2j =^3(o|e--^). (3.49) 



-N/2 

When ^ 7^ 0, and using simple properties of theta-functions [58], one may obtain instead^^ 

N/2 



hm^ J2 exp(|^(rn-Ti2)(2r-£)'j =A;^^i93(0,e--^) . (3.50) 

^ °° r=-N/2+£ ^ ^ 

As we use both results above in the ratio (3.18) for the ^-instanton partition function, we observe 
the remarkable cancelation of the elliptic/theta function contribution: the only trace of their 
existence which remains is that the result will have a different /c-dependence, depending on 
whether the instanton number is even or odd. That neither elliptic nor theta functions should be 
present in the final result is of course what one would have expected, when addressing a Stokes 
phase of a given matrix model. As such, our final result is 

Z^'^ = 1^ G, a + 1) k'-^ exp (-^) {1 + (3-51) 

where 



= Tjz- (3.52) 

2VM(a6)^/2 

In the following sections we shall test this result with great accuracy, by matching against large- 
order data. Besides the instanton action we shall give particular attention to testing the one-loop 
coefficient in the one-instanton sector (which also relates to one of the Stokes constants [7, 15]) 
which, written in terms of spectral geometry data, is very simply given by 



2V2^ (abr 



^The periodicity of the theta-function -ffs [z + mr \ q) = -ds {z \ q) implies that only the parity of £ is relevant. 
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Figure 3: Numerical data for the quartic potential. The first image shows the quartic potential V{z) — 
^z^ + j^z'^ with fi — ~1 and A = 6, while the second image displays the corresponding recursion coefficients 
r„ recursively obtained from the string equation (2.15) (in the plots with choosing N = 1000 and after 
4000 iterations in the numerical method described in [61, 62]). 

3.2 Stokes Phases and Background Independence 

In the previous subsection we used saddle-point analysis in order to explicitly find all multi- 
instanton amplitudes in a two-cuts matrix model (at least to leading order in the string coupling). 
As we have seen, the situation with a multiple number of cuts is — as long as one can evaluate all 
hyperelliptic integrals — a straightforward extension from the single-cut case [37, 7, 9]. Another 
interesting aspect of our line of work is that all these analytical results may be numerically tested 
to very high precision by making the match against large-order analysis; see, e.g., [4, 5, 6, 7, 8, 
9, 11, 12, 13, 14, 15]. As such, the obvious question to address now is whether obtaining large- 
order data for all the (generalized) multi-instanton coefficients i^^"''"^ is feasible, and perhaps 
also a simple extension from the one-cut case. In general, this is not the case and producing 
large-order data in multi-cuts situations is a much harder problem; see, e.g., [9, 14]. 

While there are several approaches to constructing large-order data, in this paper we shall 
focus solely in the orthogonal polynomial method [56] (more generally, the transseries approach 
as developed in [8, 15]). As mentioned, in general this method is in fact not applicable to multi- 
cuts configurations and what we shall discuss now is how this situation changes if we focus on a 
given Stokes phase of our system. As we also discussed in the introduction, some of the earlier 
work done in the exploration of the phase spaces of matrix models with multi-welled potentials 
was carried out in the orthogonal polynomial framework; see, e.g., [59, 60, 61, 57, 62, 63]. Such 
works were mainly based on numerical computations of the recursion coefficients, r„, appearing 
in the string equation (equation (2.15) in the case of the quartic model) and the main discovery 
concerned the appearance of multi-branch solutions at large N, as we illustrate in figure 3. 

Let us consider the case of the quartic potential V{z) = + which, when /x = —1 and 
A = 6, is depicted in the first image of figure 3. With a large N choice of = 1000 eigenvalues, 
and given the string equation for this model presented in (2.15), one may numerically iterate the 
recursion in order to compute the coefficients r„ and the result is shown in the second image of 
figure 3 (in here we have used the same numerical method as in [61, 62]). What this plot tells 
us is that, in some region of parameter space, the large N behavior of the r„ coefficients falls 
into a single branch, whereas in another region the even and odd coefficients actually split into 
alternating branches, with period two. As we shall show in the next section, this splitting of 
branches is telling us how the continuum limit should be taken in a multi-cuts Stokes phase and, 
as such, how orthogonal polynomials may be used to generate large-order results. In other words. 
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Figure 4: Numerical data for a sixth-order potential. The first image shows a sixth-order potential, 
while the second image displays the corresponding recursion coefficients r„ recursively obtained from its 
string equation (with N — 1000 and after 4000 iterations in the numerical method described in [61, 62]). 

if the recursion coefficients have a periodic large behavior, the free energy will have a well- 
defined topological expansion with exponentially suppressed instanton corrections — characteristic 
of a Stokes phase — and orthogonal polynomials may be simply used. Furthermore, notice that 
the variable n/N in the horizontal axis becomes the 't Hooft parameter in the continuum limit. 
In this case, note that the two branches merge near n/N = 1/4 which in the continuum language 
corresponds to Xt = 3/2. This critical point actually occurs when the two cuts of the quartic 
matrix model collide, and at this point the system is described in the double-scaling limit by the 
Painleve II equation. We shall have more to say about this in a later section. 

It is important to distinguish the Stokes phase, where the free energy has a "good" large 
N 't Hooft expansion, from more complicated cases which may also appear as transitions occur 
to other phases. For instance, a different behavior is shown in figure 4, obtained from the 
string equation of a sixth-order potential. We no longer find just periodic behavior, but also 
regions of quasi-periodic behavior (as shown in [37]): this quasi-periodicity is a sign of the 
theta-functions which control the recursion coefficients in this phase and which appear as one 
constructs the grand-canonical partition function of the matrix model as a sum over all choices of 
filling fractions [37]. This was recently made explicit in [41, 10], with the construction of general, 
nonperturbative, background independent partition functions for matrix models and topological 
strings in terms of theta functions. In this case, the free energy has an asymptotic large N 
behavior which is also controlled by theta functions and a naive use of orthogonal polynomials 
will not work; rather one has to use the full power of resurgent transseries. 

In summary, one may be faced with at least two different phases or backgrounds when 
addressing multi-cuts configurations: either periodic or quasi-periodic behavior of the recursion 
coefficients, corresponding to either Stokes or anti-Stokes phases. In the Stokes phase the large 
N asymptotics is essentially given by an 't Hooft topological genus expansion, while in the anti- 
Stokes phase the asymptotics is of theta-function type. These issues were addressed in [12] and 
we refer the reader to their excellent discussion (where the authors of [12] used the terminology 
of "boundary" and "interior" points to denote what we here call Stokes and anti-Stokes regions). 
In particular, an expansion around a given background is well-defined when either [12]: 

1. In a Stokes region, one will find an admissible large 't Hooft genus expansion in powers 
of with exponentially suppressed multi-instanton corrections, if 




(3.54) 
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2. In an anti-Stokes region, the free energy will display theta-function asymptotics [41, 10]. 
This expansion will be admissible if the filling fractions are real, ^ G M, and if 



The conditions of admissibility were first discussed in [17, 18], and later further addressed in 
[64, 65, 66] where they were shown to be equivalent to having the spectral curve as a Boutroux 
curve. Let us now stress that our construction in the previous subsection precisely fulfils the 
first condition above. In fact we were able to find a well-defined (exponentially suppressed) 
multi-instanton expansion, which is clear both from the general structure of (3.16) as well as 
from our final result (3.51). In this process, the Z2 symmetry plays an important role since 
it is the equality of the two instanton actions what allows us to write down a multi-instanton 
expansion for the (grand-canonical) partition function. Of course we still must make sure that 
the examples we shall address next also satisfy this condition. 

4. Large— Order Behavior of Z2— Symmetric Systems 

Our next goal is to illustrate how the multi-instanton effects we have uncovered in the previous 
section make their appearance in different examples, and how we may test them by comparing 
against large-order analysis. We shall first address the quartic matrix model in its two-cuts 
Stokes phase, as this is a particularly clean application of all our nonperturbative machinery. 
However, it is also important to have in mind that not all nonperturbative effects arise from 
what we may call i?-cycle instantons [7], i.e., instantons whose action is given by a i?-cycle 
integration of the spectral curve one-form as in figure 2. In fact, in same cases one needs to 
consider A-cycle instantons instead [11], i.e., instantons whose action arises from integrating the 
spectral curve one-form along an A-cycle and thus, because of (2.7), instantons which have an 
almost "universal" structure. As such, we shall illustrate this possibility with another example: 
the "triple" Penner matrix model which appears in the context of studying four-point correlation 
functions in the AGT set-up. Finally, notice that one of the key points that allowed us to solve 
for the nonperturbative structure of a multi-cuts configuration in the previous section was its Z2 
symmetry and, as such, this will be a required ingredient also for our following examples. 

4.1 The Two Cuts Quartic Model in the Stokes Phase 

Let us begin by addressing the quartic matrix model in its two-cuts Stokes phase. This is 
accomplished by considering the matrix model partition function (2.1) with quartic potential 



where we shall choose /U = — 1 without any loss of generality (this potential was depicted earlier, 
in figure 3). We shall first fully work out its two-cuts spectral geometry and use this data to 
obtain explicit formulae for all the nonperturbative quantities we addressed earlier in subsection 
3.1. Then, we will use orthogonal polynomials and resurgent transseries in order to, on one hand, 
readdress the results of subsection 3.1, and, on the other hand, produce large-order data that 
will be used to test and confirm our overall nonperturbative picture. 
Beginning with the spectral curve (2.4), it is simple to compute 
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Figure 5: The real part of the holomorphic effective potential for the two-cuts quartic matrix model, 
(4.4), both on the real axis (left) and on the complex plane (right), when t = 1 and A = 0.42. The brown 
areas indicate when Me Vh;eff(-z) > and the blue ones when 'ReVh-csiz) < 0. The horizontal black lines 
are precisely the cuts of the spectral curve and the black dot the pinched cycle. 

from (2.5), and the endpoints of the cuts follow from the asymptotic constraints (2.6) as 




and 



Integrating the spectral curve, the holomorphic effective potential (2.8) follows: 




(4.3) 



V^M^) = ^ I (2.2 -a^- 62) V(z2-a2)(z2_52) _ (52 _ ^2)2 / + 



The real part of this potential is shown in figure 5 where the symmetric cuts and the pinched 
cycle are very clearly identifiable. Given this result, one may immediately compute the instanton 
action, with either (3.4) or (3.38), as 



A{X,t) = Vh;efr(0) - yh;eff(a) = 1 - X " °^ ( 



(4.5) 



In its domain of validity, < At < |, this action is indeed real positive as expected. 

Similarly to what was done in the one-cut case with the quartic matrix model [7, 15], one 
may now test all our nonperturbative formulae against large-order data in a simple and explicit 
example. Of course one first needs to generate the large-order data itself and, for the present 
two-cuts scenario, the procedure will be slightly more involved than the one in [7, 15] (which, on 
what concerned the perturbative sector, was a simple extension of the pioneering work in [56]). 
Let us also stress that because this data precisely constructs the large N expansion in this phase, 
it will further confirm that it is in fact of 't Hooft type, i.e., a Stokes phase. The analysis starts 
by addressing orthogonal polynomials in this model, whose string equation (2.15) is currently 
written as 

A 
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Recall from our review in subsection 2.2 that, in the one-cut case, the recursion coefficients 
approach a single function TZ{x) with genus expansion (2.18) in its perturbative sector. This 
function satisfies a finite difference equation, (2.17), which was solved using resurgent transseries 
in [8, 15]. The key point here is that transseries solutions allow for an inclusion of all multi- 
instanton sectors, as we briefly mentioned in (2.21), going beyond the usual large N expansion. 
Furthermore, the free energy follows as (2.20). This time around, with two-cuts, as we discussed 
previously and plotted in figure 3, a numerical solution of the above recursive equation (4.6), 
approaches, in the large N limit, two distinct functions. Thus, what one now has to do is to 
generalize the aforementioned framework into a period two ansatz, as first suggested in [59, 60, 
67, 63] . As such, we shall consider 

r„ — )• V{x), n even, (4.7) 
Tn Qix), n odd. (4.8) 

In this case, the large N limit of our recursion (4.6) will split into two coupled equations 

V{x) |-1 + ^ {Q{x - gs) + r{x) + Q{x + 5,))| = X, (4.9) 

Q{x) |-1 + ^ {r{x - gs) + Q{x) + V{x + 5,))| = x, (4.10) 

and these are the equations we wish to solve via transseries methods, following the work in [15]. 
Two Parameters Transseries Solution to the String Equations 

The simplest approach to solving the above string equations, (4.9) and (4.10), is to start with a 
perturbative ansatz for both V{x) and Q{x) of the type (2.18), generalizing the work in [56], as 

+ 00 +00 

c^Y.9s'P2gix), Q{x) ^ Y.g''s'^Q2g{x). (4.11) 

9=0 g=0 

At genus zero, for instance, it is then simple to obtain 

^o(x) = 1(1-^/1-^1, (4.12) 




Qo(x) = - 1 + a/I-— , (4.13) 



where we have assumed that V Q, i.e., explicitly imposed the period-two ansatz [60, 67, 63]. 
In the domain of validity of the two-cuts phase, < Ax < |, this in fact corresponds to two 
distinct (real) functions which meet at the (critical) point Ax = |, where Po = j = Qo- 

Going beyond the perturbative large N expansion, one is first required to include all multi- 
instanton sectors via an one-parameter transseries ansatz [8, 15], 

+ 00 +00 

r{x) = ^a"p(")(x), P(")(x) ~ e-""^(^)/^^ Y,9sPt^ix), (4.14) 

n=0 g=0 
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Qix) = j;a"gW(^), qW(^) ^ ^-nAi.)/g. Y.9lQf\^)^ (4-15) 

n=0 9=0 

where we have imposed that both transseries expansions have the same structure, in particular 
that they have the same instanton action. This may a priori seem as an unnecessary assumption, 
but it is justified on two levels. On the one hand, this is required so that we may actually find non- 
trivial solutions to the string equations (4.9) and (4.10) (which are being solved "perturbatively" , 
i.e., as an expansion both in powers of the string coupling and in powers of the transseries 
parameter which corresponds to the instanton number). On the other hand, as we shall see 
later on, our large-order analysis will show that the perturbative sectors P^^\x), Q^^\x) are 
indeed governed by the same instanton action, thus "experimentally" confirming this assumption. 
Plugging these expressions back into the string equations, (4.9) and (4.10), one finds, at first order 
in instanton number and zeroth order in the string coupling, an equation for the instanton action 
as 

cosh^ iA'ix)) = 4-- (4.16) 
2Xx 

Notice that there are four sign ambiguities in this equation: two from the quadratic power and 
two from the (even) hyperbolic cosine function. For the moment we shall assume the quadratic 
sign ambiguity arises as an artifact of the period-two ansatz, and thus only address the coshz 
sign ambiguity (which is now equivalent to the one in the one-cut case [8, 15]) — and thus leave 
the complete exploration of the four sign ambiguities for future work. In this case one obtains 
for the instanton action: 



, V9 - 6Xx , / / 3 \ 

A[x) = ± — =F X arccosh I W I +27nxp + Cmt, (4.17) 

where p £ Z. We shall set both the integer ambiguity p and the integration constant cint to 
zero so that later this result will yield the Painleve II instanton action, in the corresponding 
double-scaling limit. As to the sign ambiguity, notice that choosing the upper sign makes this 
expression precisely match the instanton action as computed via spectral methods, (4.5). 

However, as shown in [15] in the one-cut case, both signs of the instanton action (4.17) are 
important when performing the fully nonperturbative resurgent transseries analysis. A similar 
situation will happen in the present two-cuts scenario, as we shall adopt the following two- 
parameters transseries ansdtze for the full nonperturbative content of the two-cuts quartic matrix 
model: 

+00 +00 +00 +00 

^(^) = ^Y1 ^1^^ P^''^"'Hx), Q{x) = ^Y1 ^1 Q^'^^'^Hx), (4.18) 

n=0 m=0 n=0 m=0 

where each p("l'")(a;) sector (and similarly for Q("l™^(x)) has an expansion of the form: 

+00 

p{n\m)^^^ ^^-{n-m)A{x)/g, ^ gfPj"l™)(x). (4.19) 

9 — i^nm 

As one plugs these expansions back into the string equations, (4.9) and (4.10), one can equate 
the terms with given powers and and find the following two coupled equations 

x6no6mo = -P^''^"'Hx)+ (4.20) 
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, n m 

A 



ni=0 mi=0 

X<5„o5n^O = (4.21) 
, n m 

ni=0 mi=0 

If one next expands these equations in powers of the string couphng, gs, this will produce — at 
each order — systems of either algebraic or (linear) differential equations which allow us to find 
the coefficients pj"''™^(x) and Q^J^^"^\x) in terms of the "earlier" ones P^^ \x) and Q^? '™ \z) 
with n' < n, m' < m and g' < g (and their derivatives). As a technical aside, let us also note that 
the many exponentials appearing in (4.20) and (4.21) via (4.19) will bring down extra powers of 
the string coupling. In fact, we shall always have in mind the following expansions: 



X +00 ^ / +00 

^ £'=0 ■ V 1=2 




(4.22) 



From here on, the extraction of the pj"'™"^ and Qg"''"^ coefficients is absolutely straightforward 
with the help of a computer, very much in line with the strategy used in [15]. Most of our explicit 
results are collected in appendix A, but for completeness we next discuss a couple of examples. 

Consider the purely perturbative sector, corresponding to n = = m, which we have also 
addressed a few paragraphs above. At order g^ it is simple to see that, once again, one finds the 
solution 



P^'^'\x) = ^"^^"^^^ = (4.23) 
A A 



Qr^(.) = ^±4^-^- (4-24) 
A A 

Here we have defined p = VS — 6Ax, as rewriting and solving most equations in terms of this 
variable will make life much easier. The remaining perturbative coefficients are recursively ob- 
tained from algebraic equations and this is generically the case for most of the [n\m) sectors (see 
the appendix A for further details and explicit expressions). 

One exception to the aforementioned straightforward algebraic procedure is when n = mil. 
In this case one finds the phenomenon of resonance, also discussed in the present context in 
[13, 15], and one needs to solve a (linear) differential equation instead. Let us illustrate this 
situation in the one-instanton sector (1|0). One finds, at order g^, 

pirn ^ 3-p ^^^^ ^^^^^^ ^(i|o) ^ ^^25) 
o 

Qi''°^ + ^cosh(A'(x))p(^l°)=0. (4.26) 

These two equations do not allow us to solve for both Pq^'''^ and Q^q^^\ but only for their ratio 

Pq^'^'*/Qq^'*^\ On the other hand, eliminating Pg^'"'' and Q^q^^\ one may instead find a differential 
equation for the instanton action — which we have solved earlier in (4.17). Proceeding to next 
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order, gl, the equations read 



P 



1 



.(i|o) , (3-p)p 
1 + 9A 
,(i|o) , (3+p)p 



sinh {A'{x)) qS,^'°^'(x) + 



sinh {A'{x)) pI^^^^^'{x) + 



3 — p 
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3+p 
6 



cosh(A'(x)) (2Q5^I°) = 0,(4.27) 

cosh(^'(x)) (2Pi(^l°^ - Po(^l°U"(x)) = 0.(4.28) 



The situation is the same as in the (1|0) sector at order g^. All we can now do is to eliminate 
the ratio P^^^'^VqI^'"^ and use our knowledge of the lower sectors — namely the relation between 
Pq^''^'* and Qq'^^\ and the result for the instanton action — in order to obtain a linear differential 
equation yielding 



These examples show a feature which is characteristic of resonance and of the n = m it 1 sectors, 
namely, that the equations we obtain at order g'^ produce differential equations whose solutions 
yield the instanton coefficients at order k — 1. At this stage the reader may object that the 
differential equations alone are not enough if one does not specify boundary conditions. In fact, 
all integration constants involved in this procedure must be fixed by using data available in the 
double-scaling limit and we shall postpone that discussion for the next section (although we have 
already used this fact in fixing the integration constants in (4.29) above). 

Other interesting features appear in the higher multi-instanton sectors, and many of these 
were first uncovered in the one-cut example studied in [15]. For example starting in the (2|1) 
sector, logarithms make their appearence into the game and they recursively propagate to the 
ensuing higher sectors. Akin to what happened in [15], these logarithms are indeed expected in 
the construction of the transseries solution and, again, we shall further discuss this issue in the 
next section, within the analysis of the Painleve II equation. Another interesting feature happens 
when n = m (and the exponential term cancels). In this case, we find that all the coefficients 
pM"') (respectively Q) with odd g vanish, and the perturbative expansion in (4.19) contains only 
powers of g^, i.e., it is an expansion in the closed string coupling. As aforementioned, further 
data is presented in appendix A, where we also find general patterns for the multi-instanton 
coefficients and relate the logarithmic sectors with the non-logarithmic ones. 

The Nonperturbative Free Energy and Large Order Analysis 

In order to test the multi-instanton results obtained in section 3, one needs to match them against 
the large-order behavior of the free energy, and this is what we shall now address. As such, we 
will derive the nonperturbative free energy of the two-cuts quartic matrix model out of the 
transseries solution to the string equations (4.9) and (4.10) we have just obtained, even though 
we will not be interested in extracting as much data. The starting point in this construction 
is expression (2.13), which yields the partition function in terms of the orthogonal-polynomial 
recursion coefficients r.„. Since in the present configuration these recursion coefficients split into 
two different branches at large A^, it is useful to first rewrite (2.13) for 2A^ eigenvalues (and thus 
with 't Hooft coupling t = 2Ngs) as 




and 




(4.29) 



2N N N 



27V-(2i-l) 

2i-i 



(4.30) 



1=1 i=l j=l 
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Notice that these equations involve derivatives of P,)^'°'(a;) and Qq^'"' (3;). 
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Similarly to what was done in (2.16), the free energy follows by taking the logarithm of the above 
expression (and normalizing against the Gaussian weight, as usual). One finds: 
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(4.31) 



It is now clear the reason why we rewrote the partition function (2.13) as (4.30) above: because of 
the even/odd split in (4.7) and (4.8), the large limit of (4.31) will precisely construct the free 
energy out of V{x) and Q{x). In the continuum limit the first sum in (4.31), which we will denote 
by the "even" sum, is essentially the same as the sum in (2.16) and thus may be computed via the 
Euler-MacLaurin formula (2.19). The second sum in (4.31), the "odd" sum, is a bit more subtle 
and requires slight modifications. In fact, from (4.24), recall that lim.x-^o Q^^^\x) ^ making 

Q^^^\x)/x ill-defined at the origin (alongside with its derivatives), but this problem is solved 
by simply considering the Gaussian contribution separately in the "odd" sector. Furthermore, 
the "odd" Euler-MacLaurin formula is now written as (following a similar analysis in [68]) 
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(4.32) 



Assembling all contributions together, our formula for the free energy finally takes a familiar 
form [56, 36, 8] 
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(4.33) 



The function iS{t,gs) comes from the Gaussian normalization in the "odd" part and is given by 
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E 

k=l 



(2iV-2f + l)log((2i- l)g,). 



(4.34) 



When computing the free energy, this expression may be first evaluated exactly and then ex- 
panded in powers of the string coupling. 

Let us note that while at the perturbative level, i.e., when n = = m, the Euler-MacLaurin 
recipe (4.33) is an efficient way to produce large-order data, the same is not valid when ad- 
dressing the (generalized) multi-instanton sectors (more on this next). In any case, using the 
expansions (4.19) when n = = m (which we have described how to compute in the paragraphs 
above, and whose data we have presented in appendix A) and inserting them into a Mathematica 
script encoding the Euler-MacLaurin expansion, we have computed the coefficients J^^*^'^^ in the 
perturbative free energy of the Z2-symmetric two-cuts quartic matrix model up to genus g = 20 
and some partial results are presented in greater detail in appendix B. 
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In order to obtain data concerning the higher instanton sectors in an effective way, and 
while remaining within the orthogonal polynomial framework, one uses a small trick due to [8]. 
Starting off with the partition function, written as either (2.13) or (4.30), it is simple to show 
that (subscripts in the partition function indicate the total number of eigenvalues considered) 

-2 = r2N+i r2N-i, (4.35) 

'^2N 

which, at the free energy level, may be written as 

T{t + 2^,) - 2J^{t) + F{t - 2gs) = log (Q {t + g^) V {t) Q{t- g^)) . (4.36) 

This expression is, in fact, a rewriting of the Euler-MacLaurin formula (4.33), but from a com- 
putational point -of-view it also makes it much easier to extract large-order data. 

We may now finally address tests of our multi-instanton formulae using large-order analysis, 
and further compute Stokes coefficients for the problem at hand. The main quantity we wish 
to focus upon is the one-instanton, one-loop coefficient Jg^'°\ At this stage, its calculation is 
simple if we are to use (4.36) above: all one has to do is to plug in two-parameters transseries 
ansdtze for all quantities and it quickly follows that, for n = 1, m = and at order (7^, one has 

4sinh^ (A'(.)) = 2 + cosh (^'(x)) . (4.37) 

If we plug in our results for the perturbative contributions, (4.23) and (4.24), for the one- 
instanton contributions, (4.29), and for the instanton action, (4.17), we finally obtain 



As we have discussed in detail in 2.3, a key point about this quantity is that it controls the leading 
large-order growth of the asymptotic perturbative expansion, as explicitly shown in (2.30). For 
completeness, let us just recall that expression in here: 

~ C l^r + ^1*' + ■ ■ ■ I ■ (4.38) 



9 ivr A^9+b I u ' 2g + b-l 

Many large-order tests may now be carried out; let us here mention a few of those following 
[7] (but, let us note, many more higher-precision tests may be carried through, as in [15], and 
these we leave for future work). One obvious test concerns the instanton action, which may be 
numerically extracted from the sequence: 

"r'-^lLr-^fi + ^ + '-V (4.40) 



The parameter b will be equal to —5/2, but that can be tested as well, e.g., using the sequence: 



6~^(25(A2af)-l)-l)+.... (4.41) 
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Figure 6: The first image depicts a test of tlic instanton action using tlie sequence ag and its first 
Ricliardson transforms, wlicn A = 0.5 and p = 1.2. Tlie iarge-order convergence towards tfie correct result 
is clear (the dashed horizontal line is the analytical prediction), with an error of 10~^% after just four 
Richardson transforms. The second image shows a test of the instanton action at fixed A = 1 but with 
varying p, after implementing four Richardson transforms. Large-order data makes up the dots, while the 
analytical prediction is given as the solid red line. Again, the match is extremely clear. 



Finally, one approach to testing the one-loop coefficient is to use the sequence: 



"^(0) r{2g + b) + • 

We should note that all sequences above have been built with free energy quantities but, of 
course, one may also perform the exact analogue large-order tests directly using the solutions to 
the string equations, V{x) and Q(x). In fact, all these quantities have their large-order behavior 
dictated by the very same instanton action^^ and, as such, we shall use either V{x) or Q{x) 
whenever possible as we have obtained far more large-order data for these quantities than for 
the free energy. We shall denote those corresponding sequences with the respective superscript. 
We also note that all these quantities have "closed string" expansions {i.e., in powers of g'^) in 
their (0|0) sectors, so the sequences above are tested for even g. 

The first natural test to do concerns the instanton action, which is shown in figure 6. Clearly, 
there is a very strong agreement between the "theoretical" prediction (be it from either saddle- 
point (4.5) or transseries (4.17) approaches) and the "numerical" data. On the left of figure 
6 we have plotted data at a particular point in moduli space^^, namely, A = 0.5 and p = 1.2, 
concerning the sequence and its first sequential Richardson extrapolations (see, e.g., [7] for 
a short discussion of Richardson transforms and their role in accelerating the convergence of a 
given sequence, within the present matrix model context). That the large-order data approaches 
the analytic prediction is very clear: after just four Richardson transforms the error is already of 
the order 10~^% at genus g = 60. On the right of figure 6 we have fixed A = 1 but vary p over 
its full range. Once again we check that the numerical data (the black dots in the figure), after 
just four Richardson transforms, is never further than 10"'^% away from the analytical prediction 
(the solid red line), thus fully validating our results. 

As we move on to testing the one-instanton, one-loop coefficient, it is important to first recall 
that the transseries framework only predicts large-order behavior up to the Stokes factors — in 



Tills was previously shown via the string equations, but we also checked it numerically to very high precision. 
^^Recall the domain of validity of the two-cuts Stokes phase, < Aa; < f , or, equlvalently, < p < 3. 
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this case up to the Stokes factor s[^\ see (4.39). However, we also have computed the same 
quantity via spectral curve analysis (3.53) (this was one of the main results in section 3.1) and, 
following [7, 8, 15], the spectral curve result should provide for the full answer. Stokes factor 
included. In this case, the calculation of S^^ Fq^^^^ in (3.53) and the calculation of F^^^'^^ in 
(4.38) combine to predict the Stokes parameter as 

. -i/J. ,4.43, 

It is quite interesting to compare the result for this "simplest" Stokes constant (at least that one 
constant which may be analytically computed from saddle-point analysis), in the present two- 
cuts configuration, with the analogue Stokes constant for the one-cut configuration in [8, 15]. 
For the quartic matrix model one thus finds: 



=-^/2 5^°) . (4.44) 

two-cuts one-cut 

With the knowledge of this Stokes constant (which we should more properly denote by sj*^^"^ 
since it refers to the free energy), we can proceed to test the relation (4.42) for the sequence 
Pg. Since besides the free energy the quantities ^^{x) and Q{x) also obey a relation similar to 
(4.42), a natural question to ask is whether the Stokes constant for these different quantities is 
the same. Indeed we find that it is the case, namely that 

This is to say that, when testing the asymptotic relation (4.42) for either f3g^^ or 

we find that the relation holds to very high accuracy with the Stokes constants being the same 
in all three cases. On the other hand, the value of b is different, with b = —1/2 for (3g^^ and 
Pg^^ and b = —5/2 for /3g"^^ (see [15] for a discussion of this point). With this knowledge, we 
have tested our instanton predictions with the sequences I3^'^\ finding that the numerical data 
has an error smaller than 10~^% at genus (7 = 60 as compared to the analytical prediction for 
S'l'^^Pg^''^^ (or Q), within most of the allowed range for A and the variable p. Note, however, that 
Pq^''''' (and also Q) diverges as one approaches p — )■ 0, making the convergence of numerical data 
to analytical prediction naturally a bit worse once we get too close to p = 0. These results are 

illustrated in figure 7. On the left of this figure we have fixed A = 0.5 and p = 1.5, and plotted 

CP) 

the sequence f3g alongside with its Richardson transforms. It is again very clear how the data 
approaches the analytical prediction (the horizontal dashed line). On the right of figure 7 we 
have fixed A = 1 and changed p over its full range, plotting the fourth Richardson transform of 
the sequence (3g (black dots) and the analytical prediction (solid red line). The agreement is, 
once again, evident. Let us mention that the very same tests may also be carried out for the 
free energy. In this case, we find an equally conclusive agreement, albeit with a smaller accuracy 
(10~'^%) as we have less large-order data available. 

At this stage one could proceed along the lines in [15] and test both multi-instanton formulae 
as well as the validity of generalized multi-instanton sectors appearing via our resurgence formu- 
lae. This would involve techniques of Borel-Pade resummation and, as such, within the context 
of the two-cuts Stokes phase of the quartic matrix model, we shall leave these precision tests 
for future work. Do notice that we shall, nonetheless, test the validity of our multi-instanton 
formulae in the double-scaling limit towards the Painleve II equation in a following section. 
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Figure 7: The first image depicts a test of tiie one-instanton, one-loop coefficient using tfie sequence 

(V) 

Pg and its first Richardson transforms, when A = 0.5 and p = 1.5. The large-order convergence towards 

the correct result is clear (the dashed horizontal line is the analytical prediction), with an error of the 

order of 10~^%. The second image shows a test of the one-instanton, one-loop coefficient at fixed A = 1 

but with varying p over its full range, after implementing just four Richardson transforms on the sequence 
(PI 

Pg . Large-order data makes up the dots, while the analytical prediction is given as the solid red line. 
As in previous cases, the agreement is extremely clear. 

4.2 The Triple Penner Potential and AGT Stokes Phenomena 

As we address nonperturbative phenomena within Stokes phases of multi-cut solutions, it is im- 
portant to note that it is not always the case that instanton effects arise from B-cycle eigenvalue 
tunneling (as discussed in [7, 9] and also as developed in section 3 of the present paper). In some 
situations, one finds systems whose instanton effects are dictated by ^-cycle eigenvalue tunneling 
instead [11, 69]. For completeness of our analysis, we shall now address an example along these 
lines. As before, we will remain within the simplified realm of two-cuts configurations, with the 
equal filling of eigenvalues ensuring Z2 symmetry of the spectral curve. 

We shall address multi-Penner matrix models. The single Penner model was first introduced 
in [70, 71, 72] and its nonperturbative effects were later addressed in [11]. Extra motivation for 
studying this system arises within the framework of the AGT conjecture [48, 73], establishing 
a relation between partition functions in 4(i = 2 superconformal quiver gauge theories and 
correlation functions in 2d conformal field theories (CFT). Within this set-up, we are particu- 
larly interested in the relation to matrix models following [74], where the quiver gauge theories 
are related to multi-Penner matrix models, and where the AGT relations follow from the the 
interconnections betweens these matrix models and CFT [75, 76]. This was further studied in 
[69], in particular addressing the three-point correlation function as a Penner matrix model cal- 
culation. The results we shall obtain below follow in this very same spirit, as they similarly 
relate to the CFT four-point correlation function with a specific, symmetric choice of insertion 
points. However, all our computations are carried through exclusively from a matrix model point 
of view, and any possible applications within the AGT context will require further examination. 

The multi-Penner potential is a sum over logarithms, as 

k 

V{z) = Y,^^^^og{z-ei). (4.46) 

i=l 

In order to obtain a Z2-symmetric potential with two wells, we shah set k = 3, 62 = I {9i + Os), 
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Figure 8: Real part of the Penner potential (4.47) for d\ = 0, = 1/2, ^3 = 1, and = ^2 = = — 1- 
fj,i = a = fj,3 and fi2 = /3. The potential now reads 



An example of such a potential is shown in figure 8 (where we plot the real part of the potential — 
the imaginary part just jumps by vr at each logarithmic singularity). The choice of parameters 
in figure 8 is precisely the case we are going to address, with = 0, ^2 = 5 and ^3 = 1. In this 
case the potential is symmetric with respect to ^; = 1/2 and not z = 0, but the motivation for 
this choice is clear: when studying four-point correlation functions on the sphere it is usual to 
make three of the insertions at z = 0, 1 and 00, with the fourth varying between and 1. Herein, 
the matrix model does not "see" the point at 00, and placing the fourth insertion at z = 1/2 
gives us the Z2 symmetry we are looking for. In the end, all that distinguishes the two cases is 
a change of variables: a rescaling and a horizontal shift. 

The saddle-point analysis we introduced in subsection 2.1 applies straightforwardly to this 
case, so we can proceed and compute the endpoints of the cuts, which we here denote by C1UC2 = 
[1/2 — 6, 1/2 — a] U [1/2 + a, 1/2 + 6]. The asymptotic behavior of the resolvent gives us three 
conditions for the endpoints, one of them being redundant. The other two are 



V{z) = alog {z - Oi) + 13 log [ z 



2(^1 + ^3) j +alog(z-03) 



(4.47) 




0, 



(4.48) 



2a + /3 + 




2t 



(4.49) 



from where we find the solutions 



a 



2 



2t {t - 2a) + /3 (/3 + 2a - 2t) - 2^t {t - 2a) {t - 13) {t - 2a - (3) 

A{l3 + 2a-2t)^ 



(4.50) 



2t (t - 2a) + /? (/3 + 2a - 2t) + 2^/t {t - 2a) {t - 13) (t - 2a - (3) 

4(/3 + 2a-2t)2 



(4.51) 
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Prom this point on, the picture is different from the one we discussed for the quartic potential. 
The main difference is that now there is no eigenvalue tunneling, in the sense that one eigenvalue 
gets removed from one of the cuts and displaced along a -B-cycle to a non-trivial saddle outside 
of that cut [7]. We can check this explicitly by looking for the zeroes of M(z) (recall that the 
eigenvalues get displaced from their cut to a pinched cycle, xq, such that M(xo) = 0). For a 
general multi-Penner potential (4.46), the moment function M{z) is given by (2.5) as 

^ (z - e,) ^j(T{6i) 

Por our particular example (4.47) we find that M{z) only has two zeroes, lying inside the cuts. 
In this case there are no non-trivial saddle points, i.e., no "hills" to place eigenvalues on top of. 
Nonperturbative tunneling effects will thus have to be distinct from our previous discussion of 
the quartic potential^^. Let us see how they arise in the following. 

One may next compute the holomorphic effective potential (2.8) and find the intricate ex- 
pression (we are using the shorthand z = z — \ /2) 



V(4a2 - 1) (462 - 1) 2a6 J ^{4a^ - 1) (462 - 1) V62 - a'^ 

X I (^8aa {2a? + 26^ - l) + ^ [c? + 6^) V(4a2 - 1) (462 - 1)^ n (</., - ^m^, m) - 

-46/3V(4a2 - 1) (462 _ i)n ((/), Vm^, H + 

+ J (6^ - a^) F((^, m) (^/3V(4a2 - 1) (462 - 1) + 8a6a) - 

-4aa (46^ - l) (n((/), m, m) + n((/-, ?i2, H) |- (4-53) 

In this expression, -F((^, m) is the incomplete elliptic integral of the first kind and n((^, n, m) the 
incomplete elliptic integral of the third kind, with n the elliptic characteristic and m = k'^ with 
k the elliptic modulus (see, e.g., [58]). Purthermore, we have introduced 



/ (6-a)(z + a) _(6 + a)2 _ (l-2a)(6 + a) _ (l + 2a)(6 + a) 

V (6 + «)(^-«)' (6-a)2' " (l + 2a)(6-a)' ' " (1 - 2a)(6 - a) ' ^ ^ ^ 

In figure 9 we show the above potential (4.53) for the choice a = —1, (3 = —1/2 and t = 1/2. 

While expression (4.53) may not be extremely insightful, it suffices to show that the holo- 
morphic effective potential is a multi-valued function. This multi-sheeted structure arises from 
the branch cuts of the square roots but, more importantly and more non-trivially, from the 
branch cuts of the elliptic functions. As we shall see in detail next, this implies that in this case 
multi-instantons are associated to eigenvalue tunneling which removes one eigenvalue from the 
endpoint of a cut and then takes it back to this cut but on a different sheet [11]. In other words, 
multi-instantons are associated to A-cycles of the spectral curve [11]. This may be shown in two 
ways. On the one hand one may analyze the branch-cut configurations of the elliptic integrals in 
(4.53) and explicitly construct the multi-sheeted structure of this function in order to study its 

^^Further note that there is no tunnehng from one cut to the other as Slj2-1 '^'^ f (•^) ~ 
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Figure 9: The real part of the holomorphic effective potential for the two-cuts "triple" Penner matrix 
model, along the real axis. In here we have chosen t = 1/2, a = —1 and /? = —1/2. 



monodromy properties. On the other hand, one may use orthogonal polynomials to first exactly 
evaluate the partition function of the model and then perform a semiclassical expansion, where 
one will be able to identify the instanton actions with j4-cycles via (2.7). For simplicity, we shall 
choose the second approach where we will find many different instanton actions reflecting the 
many different spacings between the several sheets. All these actions will be multiples of 27ri. 

As such, moving on to the orthogonal polynomial description, we first need to find the 
measure (2.11) for our "triple" Penner potential (4.47). It is simple to find 



dn(z)=e ss^^^^ — = \z\ 1/2| ssl^— 1| 

pv ; 27r ' ' ' /II I 

If we now do the very simple change of variables 

2 A = 2z - 1, dz — ;> dA = 2dz, 
the orthogonal polynomial measure becomes 



dz 
2^' 



d/i(A) = 2^+5^~Ml - A^ s« |Ar^ 



A dX 
2^' 



(4.55) 



(4.56) 



(4.57) 



The reason for doing this change of variables is because orthogonal polynomials with respect to 
this last measure, (4.57), are known. In fact, in the same way that in the single Penner potential 
(which is /c = 1 in (4.46)) we deal with the Laguerre polynomials (see, e.g., [11]) in here the 
relevant orthogonal polynomials are the generalized Gegenbauer polynomials [77, 78, 79]. Their 
precise definition is 

rl 

dX u;p,.(A) CiP'''\X) CiPf\X) = K 5nm, (4.58) 



with Wp^aiX) being the weight function 

1 



Wp,aiX) 



B{p+l/2,a + l/2) 



and 



B{x,y) 



T{x)Tiy) 



p,a > 



1 



r(x + y) 



(4.59) 
(4.60) 
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Going back to our multi-Penner measure (4.57) one immediately identifies 

2« I /S I / a 1 \ dA 
d^A) = 2^+51-^ b(i--,--^]w,_^_^{X)—. (4.61) 

In order to compute the partition function (2.13) the next step is to address the coefficients 
hn- This requires a few intermediate steps^^, for the polynomials need to be monic (i.e., the 
coefficient of the highest order term equals one). First, let us use the relation between the 
generalized Gegenbauer polynomials and the standard Jacobi polynomials Pn''^ [77] 

ciT\^) = r^^^ pif>-i/2.-m^2x' _ (4.62) 

[a + l/2)n 

ctllW = J/^±^ APrV2,.+i/2)^2A2 - 1), (4.63) 
where we used the Pochhammer symbol 

{a)m = Q (a + 1) • • • (a + m - 1) = ^ ""^ " ^^-^^^ 

i [a) 

Given this relation, one may immediately extract 

_ (p+l/2)Jp + a)Jp + a) 

- n!(a + l/2)Jp + a + 2n)' ^^'^^^ 

, (p+l/2)„(p + a)„+, jp + a) 

n\ [a + 1/2)^^-^ [p + a + 2n + l) 

These are not yet the coefficients we are looking for: as mentioned above, one must work with 
monic orthogonal polynomials and this is not the case for the Cn''^^ polynomials. But now one 
does know that Jacobi polynomials are normalized as 

^" = r(2n + ^ + z. + l) (A) ~ A + • • • , (4.67) 

which will allow us to normalize the generalized Gegenbauer polynomials. In fact, further taking 
into account the pre-factors in (4.62) and (4.63), we finally define the adequately normalized 
version of these polynomials as-'^'^ 

W- (^ + ^)^ r(2n + p + a) +•••' (4-68) 

^'^+^^ ^ {p + C7)n+l r(2n + p + a + l) 2n+lV J V J 

We now have complete information to find the correctly normalized coefficients As should 
be clear from the above analysis, they naturally split into "even" and "odd", where one finds, 

^®In these intermediate computations that follow we shall work with P = \ ^ and a — for shortness. 

Then, when addressing the partition function, we will reintroduce the original expressions. 

^'^Note that one needs to divide by 2" because J„(2x^ — 1) ~ 2"x^" H , so this cancels the 2" factor in (4.67). 
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/i2n+i = ;t-29« 9« -, -, . (4.71) 

V 9a 2gs ' 2 J \ gs 2gs 

Finally, one may compute the partition function for this model following (2.13). Due to the 
7j2 symmetry of the spectral geometry, and similarly to what happened in the quartic model in 
(4.30), the partition function naturally splits into "even" and "odd" coefficients (the ones just 
above) as one writes 

N/2-1 N/2-1 

Z= n ^2n n ^2n+l. (4.72) 

n=0 n=0 

As we shall see when addressing the calculation of the free energy, it turns out that it is useful to 
write the many products of Gamma functions which appear in the partition function as Barnes 
functions. This may be done with a reorganization of the products which appear in the expression 
above and, besides its definition T{z + 1) = zT{z), the use of the property 

r(.+ l)=2-y?0M. (4.73, 

This property is particularly useful for the Gamma functions containing half-integer factors; 
indeed the two terms with a /3 factor may now be rewritten as 

r f 2n - ^ + ll r f 2?i - ^ + 2 



r n-^ + - r n-^ + - =vr2-""-^^^^ ^ ^ (4.74) 



29s 

A similar reasoning may be applied to the terms whose numerators contain the combination 
J- + Using 7] = n — 7^ — the combination which appears in the respective denominators 
will become 

1\„/. 3^^ 



r 2r? + - r 27? + - 



r f 277 + = 2-16^,-6 r(47? + i)r(47? + 2)r(47? + 3)r(47? + 4) 
V 2^ r (27/ + 1)2 r (27? + 2)2 

(4.75) 

The final required ingredient is the definition of the Barnes function, G2{z), in terms of products 
of Gamma functions. Essentially, we shall use 

n=0 

alongside with its useful extension 

r( kn. -\- T -\- ^^ r( kn -i- T. -\- ^.^ ■ ■ ■ r( k.n -\- ^ -i- k.) = — 

G2{X+1 



n r(fc77 + X + 1) T{kn + X + 2) • • • T{kn + x + k) = '^2(fcA^ + :r + 1) ^ 

^ ^ I TO ( T" -)- 1 I 

n=0 

Assembling all different pieces together, and introducing the 't Hooft coupling t = Ngg as usual, 
the partition function of the "triple" Penner model finally follows as 

TTSs \ 



25s/ V 9s J V 25, 
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This is, of course, an exact result; it encodes both perturbative and nonperturbative con- 
tributions. In order to analyze instantons in this model we will have to understand the usual 
large N semi-classical expansion of the free energy. But this is actually a simple calculation 
given (4.78), as all one needs to know is the asymptotic expansion of the logarithm of the Barnes 
function, i.e., 

logG,(. + 1) logz - + i.log2. - ^ log. + C'(-l) + E ^^(f^ (4-79) 

As always, we are interested in the normalized free energy = F — Fq. The genus g free energies, 
J^g{t), then follow from the logarithm of the partition function (4.78), given the asymptotic 
expansion (4.79). As we are mainly interested in comparing the large-order behavior of the 
perturbative expansion against instanton data we shall only present results with genus g > 2, 
although J^Q and J-i also follow straightforwardly from this procedure. One finds: 

^^^'^ = 2g{2g-2) {'''^' " + ^''•'"^ 

+ {2^3-1 _ 1) ^2-23 ^ p2^2g _ _ p^2-2g _ _ 3^ _ ^)2-25 ^ - 2a - fif''^^) \ . 



This final result for the perturbative genus g free energies confirms that indeed the tunneling 
effects are not associated to non-trivial saddle-points as in [7] but rather to A-cycle effects as in 
[11]: in fact, its large-order growth is essentially dictated by "Bernoulli numbers growth" which 
immediately indicates that all different actions will be multiples of 27ri [11]. 

The final step we have to address is, thus, the explicit construction of the A-cycle instanton 
contributions. This again is done very much in line with the discussions in [11, 69] where each 
Barnes function in (4.78), or each Bernoulli component in (4.80), leads to a discontinuity of the 
type 

Di.c.ogG.(A'+l) = .g(M + ^)e--l-l™ (4.81) 

m=0 ^ ^ 

at the Stokes line N = \\N\. This is an exact expression for the discontinuity of the Barnes 
function and, as we discussed earlier in subsection 2.3, it immediately yields the full multi- 
instanton content of the associated free energy. In particular, the full discontinuity of the free 
energy (4.80) is given by the sum of the discontinuities of the logarithms of Barnes functions, 
and this encodes its full multi-instanton structure. We finally find 

_ 1 { { "^t ^gs\ _2li*2i 1 / 2TTt igs\ 27ritn 1 ('2'Ka l^s \ _ 27rian 

Disc-F = > < 9 6 9 6 3s e 9s -\- 

irgs IV^^ 2 \ n n"^ J 2 \ n n'' J 

'ttB igs\ 1 /27r/3 \gs\ _^2n§j± / nit - 2a) igs\ _ zeK^hMh 

+ ( 2^ ~0 2^ + 2 I ^ 
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- /3) _ igs 
n "n? 

7r(t -2a- P) 
n 

7r(2t -2a- 13) 
n 



_ 7ri(f-/3)n 

e as 



1 (2TT{t- 13) 

2 



15^ 



_ 27ri(t-^)n 
e 3a 



i{t~2a-P)r, 



+ 



n 

1 /27r(i - 2a - /3) 



(4.82) 



n 



ig^ \ 27ri(t-2a-/3)n 
^ I e 9s + 



1 /27r(2t-2a-/3) 

2 V n 



15^ 



_ 27ri{2t-2Q-/3)n 
e Us 



An alternative approach to obtaining the above result would be to first compute the Borel 
transform of the free energy (which we can do as we have an exact expression at arbitrary genus 
which, in particular, tells us that the pre-factor grows as {2g — 3)!). Then, when performing the 
inverse Borel transform out of our analytically continued result, one would extract its imaginary 
part as a sum over residues which would yield precisely the very same result as in (4.82) [11]. 

Overall, (4.82) yields twelve different types of nonperturbative effects, all with instanton 
action associated to an ^d-cycle as in (2.7), and all including only one and two-loops contributions 
around each multi-instanton sector, as further explained in [11]. The large-order behavior of 
perturbation theory will be controlled by the instanton whose absolute value of its action is 
closest to the origin in the complex Borel plane. In this example this is simple to understand 
without going into Borel analysis, by simply looking at the free energies (4.80). This expression 
is essentially a sum of numbers of the form x~'^^ and, as the genus g grows larger, the dominant 
term will be the one with the smallest x. Having said this, one may now test our nonperturbative 
discussion of this subsection by matching against large-order results. We shall do this by picking 
two distinct sets of moduli, namely, (a,/3,t) = (—1,-1,1/2) and {a,l3,t) = (-1,-1/3,2). In 
the first case the leading large-order behavior will be dictated by the dominant contribution in 
(4.82) which is the first one in that expression (of the form ~ t/2), while in the second case the 
leading large-order behavior will dictated by the corresponding dominant contribution in (4.82) 
which will now be the fourth term in that expression (of the form ~ /3/2). The instanton actions 
are A = ^ and A = respectively. We test the overall nonperturbative structure in figure 10 

by plotting, up to genus g = 100, the behavior of the sequence 1 / \fa^ (of its imaginary part, to 
be precise), alongside with the first three Richardson transforms, and for both sets of moduli. In 
a straight solid line we plot the two predictions. The error after the third Richardson transform 
is, in both cases, of the order 10~^%, fully validating our nonperturbative analysis. 



5. Asymptotics of Instantons in the Painleve II Equation 

The analysis in the previous section allowed us to check the validity of our one-instanton results, 
for the Stokes phase of two distinct multi-cut models. In particular, we have checked both the 
instanton action (3.38) and the one-loop one-instanton coefficient (3.53), predicted in section 3, 
to very high precision. But our saddle-point analysis also yields multi-instanton results, as for 
instance in (3.51), and the general structure of resurgent transseries solutions further predicts 
many, new, generalized multi-instanton sectors, as discussed in subsection 2.3. As such, we would 
now like to check all this multi-instantonic structure, and we shall do so within the context of 2d 
supergravity, or type OB string theory, by analyzing the Painleve II equation. This equation arises 
as a double-scaling limit from the two-cuts quartic matrix model we have previously analyzed, 
but is simpler to analyze from a numerical point-of-view than the full off-critical matrix model. 
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Figure 10: We plot the sequence 1/\J ag associated to ttie instanton action (see (4.40)), alongside with 
its Richardson transforms, for the choices (q;,/?,^) = (—1,-1,1/2) (left image) and (Q;,/3,t) = (-1,-1/3,2) 
(right image). The predictions for the leading asymptotics A = ti/2 and A — n/Z are given by the 
horizontal lines. In both cases the errors at genus g = 100 arc of the order 10^^%. 



5.1 Painleve II and Resurgent Transseries 

Recalling the discussions in subsections 3.2 and 4.1, it should be obvious that the two-cuts 
quartic matrix model has a natural critical point. This is cleary depicted in figure 3, which 
shows a critical point for the recursion coefRcents at At = 3/2. At this point a phase transition 
takes place, from the two-cuts phase to an unstable one-cut phase. In the double-scaling limit, 
this critical point is precisely described by the Painleve II equation. At the critical point, and 
referring to figure 5, the two cuts collide with each other, having the non-trivial saddle for the 
eigenvalue instantons at this collision point. In practice what this means is that the smaller 
endpoints of the cuts will vanish, with a — t- (and with the non-trivial saddle kept fixed at the 
origin, xq = 0). In terms of A, t and Qs, the double-scaling limit is defined as 

3 

gs^Q, A^Ac = -, t^l, (5.1) 

with the variable 

z = {l-t)gf^ (5.2) 

kept fixed in this limit. As mentioned, it is known that in this limit the matrix model describes 2d 
supergravity or type OB minimal superstrings, see, e.^., [49, 80, 51, 52, 54], and that the physics 
is encoded in the Painleve II equation. This differential equation precisely appears as we take the 
double-scaling limit in the string equations (4.9) and (4.10), as discussed in, e.g., [67, 49]. Let us 
quickly review this point, following [55], as this will also be important as we connect transseries 
solutions off and at criticality: start with the string equations (4.9) and (4.10), and introduce a 
scaling ansdtze for both P(x) and Q{x) [49] 

P(^) ^2(1- gy^ u{z) + gP v{z)') , (5.3) 

Q(x) ^2(1+ gl/-' u{z) + gf' viz)) . (5.4) 

Plugging these expressions into (an appropriate rewriting of) the string equations, it is simple to 
obtain in the double-scaling limit 

Au{z)v{z) - 2u"{z) = 0, (5.5) 
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2u^{z) - 8v{z) -2z = 0. 



(5.6) 



The first equation is readily solved for v{z) which may then be replaced in the second one. As 
such, one finally obtains a second-order differential equation for u{z), 

2u"{z)-u^{z)+zu{z) = 0. (5.7) 

The equation above is the Painleve II equation in the normalization used in, e.g., [80, 54]. In 
the present paper we are using a slightly different normalization, which follows with a simple 
rescaling of u and z as n — t- 2^/'^ u and z — t- 2^/'^ z. In this case the Painleve II equation becomes 

u"{z) - 2u^{z) + 2z u{z) = 0, (5.8) 

which, in particular, also matches the normalization used in [8]. The perturbative solution 
corresponds to an expansion around z ~ +oo where one has liport(-z) ~ y/z- 

From here on, the procedure to compute the resurgent transseries solution to the Painleve II 
equation (5.8) follows in parallel, step by step, with what was done in [15]. The one-parameter 
transseries solution to (5.8) was addressed in [8], from where we recall the following points. First, 
the perturbative solution to (5.8) yields 



-3/2 



(5.9) 



as the open string coupling. In this case, one may immediately write down an one-parameter 
transseries solution to the Painleve II equation of the form [8, 15] 



+00 +00 

U{z) ~ ^-uA/x ^nP ^ ^{n) ^ 

n=0 g=0 

where A is the instanton action and /3 a characteristic exponent. Then, plugging this expression 
back into the Painleve II equation, a solutions of this form exists if 

^ = ±3, /3=2- (5-11) 

As discussed in [13, 15] and earlier in this paper, when building nonperturbative solutions with 
transseries it is important to take into consideration all possible values of the instanton action — 
in fact, via resurgence, deep in the asymptotics of the solution one will find the need for both 
signs, and thus the need for the two-parameters transseries ansatz. As such, we shall now focus 
on the two-parameters case (but we will also recover some of the results in [8] along the way). 
Let us begin by writing the Painleve II equation in terms of a different variable 

w = x'/^ = z-^l\ (5.12) 

This is motivated by having found fi = \ above: in the two-parameters case the prefactors will 
not be of the simple form rr"^ but will depend on two integers, say n and m. As we shall see, it 
will be more convenient to include these contributions inside the perturbative expansions and, as 
such, to work directly with the variable . For simplicity of the calculation, it is also convenient 
to remove the overall factor of z^l'^ in front of the solution. This motivates us to introduce the 
new variables (with a slight, but obvious, abuse of notation) 



_ u{z) 
u[w) = — 



(5.13) 

-4/3 
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It is then a straightforward exercise to rewrite the original equation in terms of this new function 



—w^ u"{w) + —w^ u'{w) - 2u\w) -("^-2] u{w) = 0. (5.14) 
16 16 \ 4 / 

Our goal is to solve this equation with a two-parameters transseries ansatz, along the lines in 
[15], as (we remind the reader that w'^ = x is the open string coupling) 



+00 +00 

u{w,aua2) = Y.Y1 al^are-("— )^/-' <^h„)(^«). (5.15) 

n=0 m=0 



At this stage one might be tempted to assume <5(„|m)(ii;) as a power series in w but, due to 
resonance effects, this will not work (see [13, 15]): in order to obtain a solution one further needs 
to add terms multiplying powers of logu;. As such we shall use the following ansatz for the 
asymptotic expansions around generalized multi-instanton sectors 

min(n,m) min(n,m) +00 

fc=0 fc=0 3=0 

In this case, finding a two-parameters transseries solution to the Painleve II equation now trans- 
lates to determining the full list of coefficients u^"'™'^^'^'. Inserting our ansdtze (5.16) and (5.15) 
back into Painleve II (5.14), yields the recursion relation which constructs this transseries: 

n ni m mi k k~ki g g—gi 

2i: E E E EEEE 

ni=0n2=0mi=0 m2=0 fci=0 fc2=0 gi=0 32=0 



{ni|mi)[fci] „,(n2|m2)[fc2] {n-ni-n2\m.-mi~m2)[k~ki~k2] 
"91 "92 "9-91-92 



^A'in - mf + 2\ n("l™)W + '^-A{n - m){k + 1) ^"j^^t'^+'l + ^A(n - m){g - 3) u^^lf^''^ + 



+ A(^ + 2){k + 1) + ^(fc + 1)(, - 4) n("jr)['+^l + + ^(nMW_ (3^^^) 

The above recursion now allows us to see resonance explicitly. Let us consider the case where 
|n — m| = 1 and look for the leading terms in the recursion, the ■u^"'!™^^^] coefficients. The first 
term on the second line above is Qug'^^^^^\ but the sum in the first line also contains terms with 
this factor; they are: 

2^(n|m)W ^(0|0)[0] ^(0|0)[0] ^ 24010)10] ^nlHW ^OlO)^] ^ ^^OlO)^ ^(0|0)[0] ^(„|m)[fc] 

such that the leading terms in the recursion will cancel^^. As explained in greater detail in [15] 
this cancelation describes resonance in the Painleve II equation and thus the need to introduce 
the " [A;] -sectors" , which will still allow us to find a solution for the recursion in spite of the 
aforementioned cancellation. We refer the reader to [15] for further details on this phenomenon. 

Another issue which arises when solving the above recursion deals with reparametrization 
invariance of the transseries [15]: the obvious freedom to choose the parametrization of the 
transseries coefficients cji and c72 translates to a long list of free coefficients, i.e., coefficients in 
the transseries which are not fixed by the recursion. Do notice that this is not a problem, but 



""^^Recall that lirtiz^+oo Mpcrt(z) ~ ^/z, so that one has Mq"'"""' = 1. 
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rather a requirement from the transseries structure, but we refer the reader to [15] for further 
details on this phenomenon. The punch hne is that one needs to choose a prescription to fix 
these free coefficients. As it turns out, the most natural choice is to set as many free coefficients 
to zero as possible, as this will also yield the simplest final results. Following [15], we shall fix 
the reparametrization invariance by setting 

^(m+iM[o]^Q^ Vm>l, and nS"l"+'^t°l = 0, Vn > 1. (5.19) 

Having addressed the aforementioned subtleties, all one is left to do is to iterate the recursion 
in a computer. Results for the lowest sectors follow as 

^[01 , , 1 4 73 o 10657 12 13912277 

<^^,^,^Jw) = 1 w^^ w^^ , (5.20 

^ 16 512 8192 524288 ' ^ ^ 

^[01 , , 17 3 1513 . 850193 7 407117521 „ 

$ (t^) = w H H , (5.21) 

^ 96 18432 5308416 2038431744 ' ^ ^ 

1 2 41 4 5461 fi 1734407 o 925779217 m 
= -w^ + + w^^ , (5.22 

2 96 9216 1327104 254803968 ' ^ ^ 

2 291 fi 447441 ^ 886660431 14 13316458344441 

= -3w^ w^^ w^^ - ■ ■ (5.23) 

128 32768 4194304 2147483648 ^ ' 



^(210) (^) = 

^(lll)v'^ 



Let us note that, as expected, the first three lines above containing physical multi-instanton 
sectors precisely agree with the results in [8] (once we translate from our notation to theirs). 
Results concerning generalized multi-instanton sectors are new, and we present more details of 
this explicit transseries solution to the Painleve II equation in appendix C. 

We end this subsection with a few more comments on the (logarithmic) structure of the 
transseries solution and how it relates — in the double-scaling limit — to the transseries solution 
of the two-cuts quartic matrix model we have discussed in subsection 4.1 and in appendix A. The 
first thing to notice is that it is simple to determine the lowest order for which the coefficients 
^{n\m)[k] non-vanishing; let us call this number 2/3jim- The result, which can be immediately 
checked from the results above and in appendix C, is the following 



2M =n+m-2 



(5.24) 



with [★]/ denoting the integer part, and 

knm = min(n, m) - m dnm- (5.25) 

Next, and similarly to what was found for the Painleve I equation in [15], the logarithmic sectors 
turn out to be related to each other and, in particular, to the non-logarithmic sectors. In fact, 
we here find a formula very similar to the expression (5.40) in [15], which reads 

^in\m)lk] ^ i_ _ n) (5.26) 

This relation will be very useful in reducing the number of independent Stokes constants which 
enter the game; it provides relations between many of them in the same way as the analogue 
Painleve I expression was very helpful in [15]. As a final point in discussing the structure of the 
Painleve II transseries solution, let us see how to make the bridge back to the two-cuts string 
equations (4.9) and (4.10). Its two-parameters transseries solution, i.e., its coefficients p("-\^) 
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and Q^"!'") in (4.18), must agree, in the double-scaling limit, with the coefficients of our present 
solution n("l"^). That this has to be the case is clear since the Painleve II equation itself was 
derived from the aforementioned string equations via (5.3) and (5.4). But our point here is that 
this may be made explicit as we find: 

- {C VSsr- ,rV3p(n|™)[0] .-^^-^=l^nr,Sl- (5.27) 

In this expression, the "DSL" arrow simply means that we have applied the double-scaling limit 
(5.1) and (5.2) to the left-hand-side. On the right-hand-side the coefficients which appear are 
the ones associated to the original variables, i.e., where we have inverted the redefinitions of u 
and z we did before. There is a similar expression envolving the Q coefficients which relates to 
the one above by a simple change of sign, as can be seen in (5.4). Finally, the constant C is given 
by 

C = (5.28) 



A 

5.2 The Resurgence of Multi Instantons and Stokes Coefficients 

We shall now turn to the resurgence of the generalized multi-instanton (n|m)[A;] sectors. The 
resurgence formulae we have discussed in subsection 2.3 will verify the validity of these multi- 
instanton sectors, and they will further allow — upon consistency — to extract many unknown 
Stokes constants. We shall only focus on effects at exponential order and our analysis will be 
less detailed than the one in [15] where, using more refined techniques, it was possible to "dig" 
deep in the asymptotics and study effects at orders 2~^, 3~^, et cetera. Nonetheless, our results 
will fully validate the two-parameters multi-instantonic structure of the Painleve II solutions. 

Let us begin by addressing the Stokes constant On what concerns large-order behavior, 
this constant appears in the perturbative (0|0) sector and we may use the large-order expression 
(2.30) to write in the present case 

,,(0|0)[0] _ Sf^ r (2g - I) (i|o)[o] ,fe r(2ff-/i- ^) 

- ^2,-1 r(2<7-i) + ^ ^ ^ ^ 

Given this expression, it is immediate to construct the sequence 

ITtA^^-^ (0|0)[0] 



r(2,-i) 



<J"^'"' (5.30) 



which is asymptotic to 5^'^^ . Taking its Richardson extrapolation, it follows an extremely precise 
check on the well-known result (see, e.g., [54, 8]), where we found a match of the first 30 decimal 
places after = 20 Richardson transforms 

S\°^ = ^ = -0.3989422804014327... i. (5.31) 

V2^ 



There is a simple relation between the above Stokes constant at criticality, and the corresponding 
Stokes constant off-criticality, (4.43), which is similar to the relation between the corresponding 
Stokes constants in [15] — i.e., Stokes constant for Painleve I and for the one-cut quartic matrix 
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model. Namely, we find"*^^ 

^(0) 

(0) _ ^ QMM 



(5.32) 



PII C 

where the constant C was defined above. Naturally, this expression is simply encoding the 
double-scaling limit at the level of Stokes constants (see [15] for other Stokes constants). 

As we move forward there is one point to have in mind: except for a limited set of empirical 
relations they satisfy among themselves — which we shall discuss in the following — there are no 
further analytical predictions for all other Stokes constants. As such, we need to compute them at 
the same time we test resurgence in an independent fashion. This is done in two steps [13, 15]: we 
choose one resurgent formula and validate it via some resurgent relations; then we use different 
resurgent relations in this formula to numerically compute new Stokes constants. As one iterates 
this procedure towards several Stokes constants and several multi-instanton sectors, consistency 
independently double-checks both the Stokes constants and the resurgence of instantons. 

In this spirit, let us move on to the multi-instanton sectors and address the Stokes constant 
S_l which appears in the (2|0) sector. If we apply our large-order formula for multi-instanton 
sectors, (2.31), with two physical instantons, n = 2 and m = 0, and focus only on the leading 
contributions to the asymptotics, = 1, we arrive at 



u. 



(2|o)[o] _ 3^}°^ ^{a-D r (g - /i - I) (3|o)[o] 

2^+2 - 2^i ^5-1 r(g-i) ''2/.+3 ^ + i5-33j 

I ^-1 ^ (g ~ 5) V" ^js-h-^) (2|1)[0] .Nh , 

+ 2.i (_^).-| A. r(g-i) ^--3 (-^) + 

+ 4vri (_^).+ | A. r(^+i) -2H^i i ^) h+^) log (A) 



ITT 



A novel feature of this case is that, adding to the familiar g! large-order growth, the digamma 
function further produces effects which grow as g!logg and which will in fact be the dominant 
effects. One procedure to extract and confirm the new Stokes coefficients associated to this 
expression, via Richardson transforms and when in the presence of log g factors, was introduced 
in [13] within the context of the Painleve I equation and further extended in [15]. Let us see how 
to address this issue. We move a factor of 

2vri ^. ^y. (5.34) 
r (9 + 2) 

to the left-hand-side of the above equation, and expand its right-hand-side in powers of 1/g. In 
this way, one obtains a sequence with the following asymptotic behavior: 

+00 ^ +00 

Ag^ Bglogg + Cg, wheie ^<7-J^4, ^^ff-Z^S' ^^'^^^ 

k=0 ^ k=0 ^ 



^^FoT shortness we will avoid the labels referring to either "Painleve H" or "Quartic Matrix Model" throughout 
the rest of the paper. All constants discussed from here onwards refer to the critical (double-scaled) model. 
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To extract the leading coefficient, b^, we may construct a new sequence, 



Ag=g {Ag+l - Ag) , (5.36) 



which behaves as 



+00 T +00 ~ 

Ag^Bglogg + Cg, whcrc 5^-^^' ^9-^0 + ^^, (5.37) 

k=l ^ k=l ^ 

thus isolating the coefficient we are looking for. In fact, should we apply a couple of Richardson 
transform (at least two), we remove the subleading tails in l/g'^ and in log g/g^ and immediately 
obtain bo numerically. Similarly, if we want to find 6i, we can define 

=g[Ag- bo) , (5.38) 

and now apply the Richardson transforms to the sequence A^g^ = g ^^^+1 — ^g^^^ iii order to 
extract the coefficient —bi. As we move on to the extraction of the Cj coefficients, the procedure 
is more or less straightforward. For instance, if we subtract the leading logarithm to the left- 
hand-side of the original sequence, the new sequence^*^ 

Pg = Ag- bo log g (5.39) 

will now yield cq. Along the same lines, Pg^^ = Ag — {bo + bi/g) logg allows us to extract — ci, 
and so on. Applying all this in our present context we now have to consider the sequence 



21 



^. = 2vri-^ng^f, (5.40) 

where we should notice that, due to the factors of (—1)^ in (5.33), we need to look separately 
at the sequences for g odd and for g even. For simplicity, we shall only discuss the even case, 
but the odd one is completely analogous. If we now use the sequence (5.36) to compute A2g, we 
expect it to converge towards the leading coefficient multiplying logg in the resurgent relation 
(5.33). What is this number? Using the value of u^^'^^'^^ = —8 (simply obtained for instance 
via (5.26)) and using the fact^^ that S^l = —iS^\ if the resurgent formulae hold in the present 
context then this number should be equal to the analytical value 



iS 



(0) 

-1 .,{2|1)[1] 



uY^^>^'^ = 1.59576.... (5.41) 

Let us then turn to the sequence and analyze it. This is shown on the first image of figure 
11, where we plot the original sequence and some of its Richardson transforms. After = 20 
Richardson transforms we find the numerical value of 1.59573... which differs from the prediction 
above by less than 0.01%, thus fully validating our resurgent multi-instanton structure. 



^"Notice that we can subtract further logarithmic terms in order to accelerate the convergence. 
■^^A trivial word on notation: A is the instanton action, Ag the sequence we are addressing. 
^^At this precise moment this only adds numerical evidence to the fact that S^l = -15*1°^ But, as we shall see 
in the following, we can actually show that this relation is true, so we may as well use it in the following. 
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Figure 11: The left image shows tlie sequence A^g built from (5.40) (blue) alongside with its fifth (green) 
and twentieth (red) Richardson transforms. This can be shown to quickly converge towards our prediction 
(5.41) with errors ~ 0.01%. The right image shows the sequence Pg in (5.39) (blue) alongside its fifth 
(green) and twentieth (red) Richardson transforms. This quickly converges towards our prediction (5.43). 



Taking our analysis one step further, we may now extract a new Stokes constant by looking at 
the leading non-logarithmic term, which may be computed using the sequence (5.39). According 
to the same large-order resurgent relation, (5.33), this term should be 

I'^-l (2|1)[1] A \ ■ q{2) (1|0)[0] 
— u\ [^ogA + ivr) — id1{ u\ , (5-42) 

(2) 

where we do not know the value of the Stokes constant S_(. However, we may find it by analyzing 
the sequence (5.39), as shown in the second image of figure 11: after = 20 Richardson iterations 
we extract the numerical prediction 

S^y^ = -5.3455144... - 5.013256493... i. (5.43) 

Before moving on with further Stokes constants, let us make a remark concerning the new 
Stokes constant we have just computed, (5.43): it is a complex number, with both real and 
imaginary contributions. But, as explained in detail in [15], there are many relations between 
the Stokes constants and a large number of these depend on each other (although it is unclear 
how many truly independent Stokes constants exist). Some of these relations may be derived 
from the general structure of the string genus expansion, and are thus model-independent; while 
others were found "experimentally" , and will thus depend upon which equation is under analysis 
(but see [15] for more details on both these points). In particular, all Stokes constants of the 

(n) ~('n) 

form and with i > are purely imaginary. We will thus only list this type of Stokes 
constants. We shall discuss how these relations arise when we discuss the (1|1) sector below; for 
the moment let us just mention that, for (5.43) above, the relation which involves out of 

5f) = -i^i^J +47ri5j°^ = 5.3455144... i. (5.44) 

Having successfully addressed a two-instantons sector, let us next address a sector involving 
generalized instantons. In this case, the simplest choice is to study a generalized "closed string" 
sector; the example where we have n = 1 = m and A; = 0. By "closed string" we mean that 
sectors of the type (n|n) are expected to have an asymptotic expansion in powers of the closed 
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string coupling ~ rather than in powers of the open string couphng gs 
seen in (5.23) [15]. In this case the (1|1) sector has no logarithmic contributions and the relevant 
large-order relation is [15] 

27riA^^+5 (i|i)[o] (1) :^r{g-h + l) (i|o)[o] ^ , ocW 4 ^ ~ ^ ~ 5) „,{2|i)[o 



ZTTIA" 2 (i|^)[o] ^ n) ^ 1 yg-n^^) (i|o)[o] ^(0) , 1 " " ^ (2|1)[0] 

2c(0) V ^ ~ ^ + 2) ,/2|l)[l] Ahf^ ( _ 



h=0 

iff) :^ r{g-h + l) (o|i)[o] .^fe 2iffj . r (g - - ^) (i|2)[o] ^ 
2- r fo + i-^J +731)7^2^ 



+ 



As we have just mentioned, the (1|1)[0] sector will have a standard, topological perturbative 
expansion, which implies that all terms above with odd g will have to vanish [15]. In other 
words, imposing ^^2(23+1)4-2 ~ ^'"^ result in a tower of relations between the Stokes constants 
appearing on the right -hand-side of (5.45), as this expression needs to vanish order by order in 
both powers of \jg^ and logg/g^. For example, expanding the digamma functions we find that 
imposing that the term proportional to log g vanishes will imply the condition 

5W-iff) = ^ S^'} = -isf\ (5.46) 

which we had already put forward and checked numerically — now being "theoretically" justified. 
On the other hand, the term at order 0(1) yields a relation involving two unknown constants, 

S[^^ + i s'l} + Svri = 0, (5.47) 

where we have used (5.26) to relate Continuing along these lines and looking 

at further required cancellations, one may use this procedure in order to extract similar relations 
between further Stokes constants, such as (5.44) which we have discussed above. Our goal now 
is to apply the same reasoning as used within the (2)0) sector in order to compute this new 
Stokes constant, s\^^ (and, along the way, S^l as well). This is very similar to what we have 
done before, with the slight difference that now only the sequences for even g are relevant. Once 
again the term proportional to log g offers just a consistency check on the resurgent structure of 
the transseries solution and on (already) known Stokes constants, and we show in figure 12 that 
this is indeed working perfectly: the relevant sequence, after Richardson extrapolation, converges 
towards the correct number, —5^*^^ u^^'^^'^^ + iS^l u^^'^^^^\ with an error smaller than 0.001%. 

The new constants we are after appear at order 0(1), without logarithmic contributions. 
After using the relevant sequence, (5.39), we find, as we show in figure 12, a fast convergence 
towards the number ^ = —19.54576... i that resurgence sets to 

e = S['^ Wl^l log^ + S['^ n(^l°)M -iff) n(°WM - iff°) nfl^)™ (log^ + ivr) . (5.48) 

Using this result together with the previous relation, (5.47), we find s[^^ (which is, as expected, 
a purely imaginary number) 

S-W = -10.6910288... i. (5.49) 
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Figure 12: The left image shows the sequence which tests the leading, log 5, coefficient of the large-order 
relation (5.45) (blue), alongside with its fifth (green) and twentieth (red) Richardson transforms. This 
sequence quickly converges towards the expected limit — sj^-* m^^'^-"^' + i S^_^lu^™ with an error smaller 
than 0.00f%. The right image shows the sequence which tests the leading, order 0(1), term in the large- 
order relation (5.45) (blue), alongside its fifth (green) and twentieth (red) Richardson transforms. Again, 
this quickly leads to our prediction (5.48). 
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-0.39894228... i 
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-10.6910288... i 
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5.3455144... i 
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Table 1: The independent Stokes constants we have calculated. The third column gives the number of 
decimal places to which the answer is explicitly computed, while the fourth column shows the instanton 
sector where each constant appears for the first time. All constants we address first appear at order 



The (independent) Stokes constants we have computed are summarized in table 1. It is 
interesting to notice that a further "experimental" relation s[^^ = — ^ s[^^ is (apparently) true 
in this case. The exact same relation was also found in [15], in the context of the Painleve I 
equation, alongside with some other extra relations, all of them emerging from purely numerical 
relations. We expect that by examining further data in the present Painleve II context also many 
similar relations will be found. However, at this stage, we have no first principles explanation 
for these extra relations: determining the minimal set of independent Stokes constants is a very 
interesting open problem for future research. 

5.3 The Nonperturbative Free Energy of Type OB String Theory 

The final point we wish to address is the construction of the nonperturbative free energy for 2d 
supergravity or Id type OB string theory. In fact, using the results of our transseries analysis of 
the Painleve II equation, we may now build its associated double-scaled free energy. This free 
energy is obtained from the solution of the Painleve II equation via [49, 80, 51, 52] 

F:i,{z) = -\u{zf. (5.50) 

For convenience, from this point on we shall drop the double-scaled label, but we will always be 
talking about the free energy at the critical point. The first thing to notice is that there is now a 
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fundamental difference with respect to the Painleve I case, studied in [15]: the relation between 
the (twice derived) free energy and the solution of the differential equation is no longer linear. 
Nonetheless, the right-hand-side of (5.50) still has a transseries expansion 



-. +00 +00 

4^W' = E E <^2"^e-("— )^^''Vh„)(z), (5.51) 



n=0 m=0 

but where now one has 

n m +00 ~(n|m) 

'P{n\m){z) = E E '^{n'\m'){z) (n-n'\m-m'){z) ^ E ^ M • (^-52) 
n'=Om'=0 9=0 ^ * 

Relating this expression to the free energy now just requires a double-integration, as follows from 
(5.50). Let us begin by looking at the perturbative sector, where we bring back the y/z factor 
we had pulled out in (5.13). In this case, the integration leads to 

..0,0,,., . 4 11 -4-'^.^,.^..^-. (5.53) 

As a check on this result, notice that if we apply the double-scaling limit, (5.1) and (5.2), to 
the quartic matrix model free energies, Fg{t), which we have computed via the Euler-Maclaurin 
formula in (4.33) (the first few of which are presented in appendix B), and if we further implement 
the rescalings and z — )• 2^/'^ z associated to our choice of normalization, then the 

answer one obtains precisely matches the above result. 

Having understood how to construct the free energy in the perturbative sector, one may 
move on towards multi-instanton sectors. Beginning with the one-instanton sector arising from 
the product <I>(o|o)(-2) ^(i\o){z), the first coefficient to compute is simply given by 

_1^^4i|0)[0] jjdzz'/^e-^^'^' = -lainfl^^M.-^'/^e-^^^^^ + ••• . (5.54) 

In the expression above we have kept only the leading term and we have explicitly displayed the 
coefficient u[^''^^^'^'. Recall that when solving the Painleve II equation we chose to set u^^''^^^'^' = 1, 
and recall that this freedom in choosing the normalization was a consequence of a reparametriza- 
tion invariance of the double-transseries solution [15]. One now needs to readdress this point 
in order to properly fix the free energy transseries. As shown in [15], rescaling the transseries 
parameters as cri = ci ai and o"2 = C2 ^2 makes the following quantities scale accordingly 

^{n|m) = q"c2™$(n|m), (5-55) 

5f =ci-'=4-'=-^5f , (5.56) 

^(k) 

5f =ci+^-^4-'=5, . (5.57) 
The convenient scaling to do, when dealing with the free energy, is 

cTi,2 = 5r<2- (5-58) 

In fact, this immediately implies that the leading coefficient of the one-instanton free energy is 

^(i|o)[o]^_l5(o)^(i|o)[o]^ i ^553^ 
° 8 ^ ^ 8^2^ 
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and thus the free energy Stokes constants is very simply 

= 1. (5.60) 

This convenient normahzation may further be double-checked by using the large-order connec- 
tion (4.39) 

pimm ^ ^ ~ i) ^(i|o)[o] g-^^ 

After just a few Richardson transforms we find that (5.60) is indeed consistent. It is important 
to remark that physical quantities should not depend on normalization choices, so that only 
combinations which are left invariant by the above rescalings are physical. In this particular 
case, the physical quantity is 

(5.62) 

Had we chosen to have = ai, then we would have found 5^*^^^ = sf^\ but the combination 
above would not have changed. A longer discussion on normalizations may be found in [15]. 

We are now ready to proceed and explicitly compute generalized (n|m) multi-instanton 
sectors in the free energy of 2d supergravity or Id type OB string theory. From the point-of-view 
of the double-integration, the only complicated sectors are the ones with logarithms. In fact, 
when n = m the procedure is immediate and a straightforward generalization of what we did 
for the perturbative (0|0) sector in (5.53). As such, and always having in mind that we are now 
dealing with the function u{z)'^, in (5.51), we have for general n / ?n, 

e-("-'")^/'"' V\1\^){W) ^ a'^a^ ^-{n-m)AI^^ ^ ~(n|rn)[0] ^23+2/3il_ (5 ^3) 

9=0 "™ 

It can be shown — and easily checked — that the relation (5.26) connecting logarithmic (n + A;|m + 
A:)[fc] to non-logarithmic (n|m)[0] sectors still holds in the precise same form for ^'^'^ 
components. In this case, it is convenient to assemble together all sectors which are related to 
the (n|m)[0] sector; due to the aforementioned relation each of these is of the form (transseries 
parameters and logarithmic factor included) 

i (8 (m - n) CT1CT2 logu;)'^ a^a^ ^-[n-m)A/^^ V^[°i„)(^). (5-64) 

Finally summing over A;, one finds 

In order to do the double-integration, let us first move back to the z variable so that the {n\m) 
contribution becomes 

+ 00 [Q] 

_(._^)A.3/2 y~in\ram ^_5.±5^ + 5(:...^ _ 
9=0 

The double-integration may now be carried through using that 

O +00 _ 

dz z'^ e-^^^''' = ^ e-'^^'^' Yl • {-^Az^^') " > (5-6T) 

m=l 
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where the coefficients am{q) are given by 

Notice that the amio) coefficients are polynomials in q of degree m — 1. Further, given the 
integrand in (5.66), the variable q is actually linear in aia2 and, as such, the am{q) coefficients 
will be polynomials of degree m — 1 in (Ti(T2. This effectively means that the double-integration 
of the (n|m) sector of u{z)'^ contributes not only to the (n|m) sector of the free energy, but to 
all other (n + r|m + r) sectors as well (with r > 0). 

We are now essentially done. Using a computer, we can apply the integral (5.67) systemati- 
cally and find that the free energy has the structure 

+ 00 +00 

n=0 771=0 

where the "coefficients" F^"'^'^\z) will be asymptotic expansions in powers of (both integer 

and half-integer, and also containing the occasional logarithm). The first few sectors of the 
critical free energy are the following 



i?(o 


°)(z) 


1 3 

= -24^ 


1 1 3 63 c 
32^°^"+ 512" + 4096^ 


(5.70) 




°)(z) 


1 a 

= 


65 9 19273 15 13647905 21 
"^ 768^ * 147456 ' ' 42467328^ 


(5.71) 




1)(.) 


4 3 


25 3 6323 9 5015413 15 

Z 2 4- z 2 -|- z 2 4- . . . 

96 24576 3145728 


(5.72) 




°)(.) 


1 


3 59 _^_3 9745 _^_9 3335669 
' 1536^ UlAbQ'^ " ' 21233664" 


(5.73) 


i?(2 


i)(z) 


9 

= -16^" 


9 737 _i5 398375 „2i 142017823 „27 

4 -U Z 4 — Z 4 -1- Z 4 — • • • 

512 98304 9437184 


(5.74) 


i?{2 


2)(z) 


= -3 log 


111 54507 15245711 „ 

z H z~-^ -\ z~^ H H , 

64 8192 196608 


(5.75) 


i?(3 


°)(.) 


1 

~ ~96^ 


9 59 15 7645 21 1836031 27 

4 -| z 4 z 4 H z 4 — • • • 

3072 196608 18874368 


(5.76) 


i?{3 


i)(z) 


17 

= "64^" 


1211 9 655883 . 161783969 15 

H 2; 2 z " H z 2 

1536 294912 21233664 


(5.77) 


i?(3 


2)(.) 


17 9 

= Y'~' 


2267 15 3488915 21 251878099 27 
384 ' ' 147456 ' 2654208 ^ * + 


(5.78) 


i?(3 


3)(z) 


17 3 

= — z 2 
3 


35675 9 11452163 15 157674856009 21 

H z 2 -1 z 2 -1 z 2 + . . . 

2304 81920 58720256 


(5.79) 



In the list above we presented the sectors (n|m) with n > m. The coefficients with n < m differ 
at most by signs, obeying the rule 

pMn) ^ ^_iy+[n/2]j pin\m) ^ n > 771. (5.80) 

The starting powers in the free energy coefficients ^("1*") can be easily related to the starting 
powers l^nm of the Painleve II coefficients u^"!™-), for instance by looking at (5.67). At the end of 
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the day we find 



(5.81) 
(5.82) 



where /3nm was defined above in (5.24). For completeness, we also recall that the logarithmic 
sectors are "hidden" inside the term 



(m-n)af<xf = exp ( ^ (m - n) af erf logz 



(5.83) 



As a final point, we should also comment on the Stokes constants for the free energy. Since the 
Stokes constants for u{z) and u{z)'^ are the same, the Stokes constants for the free energy are 
related to those of Painleve II via the rescalings described above, (5.56) and (5.57), and so 



2k+e-2 



(k) 



«2 t qiO)^^^-^-^ g{k) _ 



ib.M) 
(5.85) 



The extra £ appearing above comes from taking two derivatives on the factor exp 

On what concerns the independent Stokes constants we computed in section 5.2, the respective 

values for the independent free energy Stokes constants are presented in table 2. 







Precision 




1.0000000000... 


oo 




-4.26510341... i 


8 


gi2)F 


2.13255170... i 


8 



Table 2: The independent Stokes constants for the free energy of 2d supergravity or Id type OB string 
theory. They are related to the Stokes constants of the Painleve II equation via (5.84) and (5.85). 



6. Conclusions and Outlook 

In this paper we have continued our analysis of the nonperturbative structure of the large limit, 
and of string theory, along the lines in [15]. We have generalized many results [7, 8, 9, 11, 13, 14] to 
the multi-cuts realm, with emphasis on proceeding with the nonperturbative study of the quartic 
matrix model initiated in [7, 8, 15], this time around in its two-cuts phase with its Painleve 
II double-scaling limit. Our results support the need for resurgent analysis and transseries, 
in particular the need for two-parameters transseries solutions including new nonperturbative 
sectors. As in previous work, the question remains to explain, semi-classically, what these new 
sectors are: for example, in the Painleve II context, while the physical multi-instanton sectors 
correspond to ZZ-branes [81, 82, 8] there is no similar understanding of the generalized sectors. 
Partial discussions may be found in [13, 14, 15, 83] but no conclusive answer as yet been reached. 
This question is also related to finding a first-principles calculation of the many "experimental" 
Stokes constants we have found: for one of these constants, in the Painleve II framework. 
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there are many analytical methods which determine it, see, e.g., [53, 47, 54, 8, 84], but for all 
others finding one such analytical method is still an open problem. This is probably related to 
first determining how many truly independent Stokes constants there are for each problem, and 
further explaining the empirical relations we have found among them. 

As one looks towards future research, some natural generalizations of our present work 
quickly come to mind. For example, one natural extension would be to "dig" deeper into the 
asymptotics of our examples, putting the full resurgent formulae on even stronger grounds. 
Recall that in [15], both for the one-cut quartic model and for its Painleve I double-scaling 
limit, techniques of Borel-Pade resummation were used in order to analyze contributions to 
the large-order behavior arising at exponentially suppressed orders of 2~^, 3~^, and so on. It 
would be very interesting to extend those results within the present examples of the two-cuts 
quartic model and its Painleve II double-scaling limit. Another interesting line of work would 
be to further explore the connection to the AGT framework, along the lines in [74, 69]. For 
example, one could address the calculation of (symmetric) higher-point correlation functions, or, 
in a different line, compare our present nonperturbative construction with other nonperturbative 
completions suggested within the AGT set-up as in [85]. Yet another interesting line of work 
would be to address extensions of our Painleve II results towards its deformations which arise 
within the type OB minimal superstring context, when turning on RR flux or in the presence of 
charged ZZ-branes [50, 52]. This flux is controlled by a parameter, q, and the equation which 
describes the minimal string set-up is now a deformation of Painleve II, namely 



This equation is certainly addressable within our framework and it would be very interesting to 
fully carry out its resurgent transseries analysis, extending our Painleve II results. 

In order to be fully explicit when addressing multi-cuts Stokes phases, we have focused 
on the two-cuts case where the Stokes phase is essentially related to the Z2 symmetry of the 
spectral curve configuration. But one may extend this calculation for an arbitrary number 
of cuts, k, as long as one keeps the corresponding spectral geometry configuration having a 
Zfc symmetry on its cuts — this symmetry will ensure that, although generically dealing with 
hyperelliptic configurations, at the end of the day all calculations reduce to elliptic integrals 
(very much along the same lines as it occurred for us in subsection 3.1). Afterwards, and still 
following our own guidelines from the Z2 case, a proper treatment of the sum over instantons will 
further ensure that these elliptic functions will cancel in the end, thus producing adequate results 
for a Stokes phase. Setting up such symmetric spectral configurations is very simple, as it is to 
compute their corresponding instanton actions. The multi-instanton analysis should then follow 
with some extra work. Another interesting point of this example is its own double-scaling limit 
[67, 86] which seems to lead to new integrable hierarchies. For a Z^-symmetric configuration the 
string equations get more complicated, but are certainly solvable within our framework. In this 
way, it should be possible to say a lot about the nonperturbative structure of their corresponding 
solutions and, thus, about the general structure of these new integrable hierarchies. 

The two-parameters transseries solution we have obtained for the Painleve II equation is, in 
principle, its full nonperturbative solution. How may we understand the information encoded in 
this solution? When addressing 2d supergravity, or type OB string theory, we are looking for a real 
solution to this equation, (5.8), for all z € M. Recall from, e.g., [8] that this is naturally associated 
to the two phases of the Painleve II solutions: the weak-coupling phase, when z — t- +00, and the 
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strong-coupling phase, when z 



—oo, where one finds the asymptotic behaviors [87, 53] 




z — )• +00 



(6.2) 



u{z) 



1 




l/4g-2^(-,)3/2 



(6.3) 



21/12 ^ 



z 



—00. 



There is in fact one such real solution, interpolating between the above weak and strong couplings, 
the Hastings-McLeod solution [87]. Notice that there are many solutions to the Painleve II 
equation: for example, in [53] one finds a large class of global purely imaginary solutions, alongside 
another large class of global real solutions of which the Hastings-McLeod solution is a particular 
(singular limit) case — and the asymptotics of all these solutions are well known [53]. In particular, 
all these solutions should be encoded in our transseries solution, but in here we shall only discuss 
the Hastings-McLeod solution which was also addressed in [8]. This in itself is already a non- 
trivial problem: clearly, the "instanton action" in (6.3) is different from the Painleve II instanton 
actions appearing in its two-parameters transseries solution. The natural question that follows 
is: how is the Hastings-McLeod solution encoded in our two-parameters transseries solution, 
and how can it provide for both types of weak and strong coupling behaviors, displayed above? 
In particular, how may A and —A of Painleve II "conspire" to yield the extra factor? This 
question was partially addressed in [8], in the context of an one-parameter transseries solution. 
In there, it was shown that — upon Borel resummation — one may perform a median resummation 
of the transseries along the Stokes line in the positive real axis to yield a real solution of the 
Painleve II equation (see the final discussion in [15] as well), i.e., 



where 5+ denotes a left Borel resummation along the positive real axis (see, e.g., [15]). Once this 
is done, the Hastings-McLeod solution is that particular real solution which has cr = in the 
expression above [8]. In particular, this median resummation of the one-parameter transseries 
reproduces the Hastings-McLeod content for z G M"*". But one question remained open: what 
happens along the negative real axis instead? To answer this question one needs the full two- 
parameters transseries solution we have constructed in this paper, but yet this is not the full 
story: constructing a median resummation along the negative real axis, where one now finds an 
infinite number of highly non-trivial Stokes constants, is a much harder problem, and moving 
from positive z to negative z will also entail crossing Stokes lines. These crossings will make 
Stokes constants jump, not only as overall factors but also inside exponential terms due to the 
logarithmic sectors as we discussed in subsection 5.3. As such, it would be a very interesting 
project to make this strong/weak coupling interpolation completely explicit, within the resurgent 
transseries framework. We hope to return to some of these ideas in the future. 
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A. The Two-Cuts Quartic Matrix Model: Structural Data 



In this appendix we present some explicit results concerning the two-parameters transseries 
solution to the two-cuts quartic matrix model. Let us recall that in subsection 4.1 we have 
solved the string equations of this model, (4.20) and (4.21), by introducing the ansatz 

+ 00 +00 

n=0 m=0 

with 

+00 

p(nM^^) ^g-{n-m)A{x)/5s ^ 5fpH™)(^)^ (A.2) 

Q — i^nm 

and similarly for Q{x). In the table below we show the maximum order (in g) in the string 
coupling to which we have recursively computed the above nonperturbative coefficients: 



\ n 
m\ 





1 


2 


3 


4 





60 


10 


10 


10 


5 


1 




10 


5 


5 


5 



Table 3: Values for the highest g for which we have calculated pj"'™-* and Qg"'™''. 

Do note that, since the sums in (A.2) have a "starting genus" which is (in general) Pnm = 
— min(m, n) < 0, the actual number of coefficients that we have computed is bigger than the 
numbers displayed in table 3. It is also worth pointing out that in the cases where n = m 
the asymptotic expansions contain only even powers of gs, which implies half of the indicated 
coefficients vanish. Finally, the sectors (?7T,|n) and {n\m) are trivially related via (similar for 
Q{x)) 

pHrn)(^) = (-i)9pH")(^). (A.3) 

Let us begin by presenting explicit results for the first few coefficients in the perturbative 
sector^^ 

Po'°'°' = i(3-p), Qr = i(3 + P). (A.4) 



2p5 ' 2p5 

,(0|0) _ ^3 1915812 - 314928p - 181521^2 + 1871 V + 1944/ 



p(U|U) _ x3 J^^J-'Jui^ j.u-L<j^i/y -r j^uij^x^ -r i^/n^ 



(o|o) _ ,3 -1915812 - 314928p + 181521^^ + 18711^^ - 1944/ 
Q2 - A . (A.7) 

Proceeding with the multi-instanton sectors (and just explicitly showing results for V{x) from 
now on), the first few coefficients in the (1|0) sector are 



pr = -f-^. (A. 



^^Recall from the main body of the text that we are using the variable p = y/9 — 6\x. 
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8p7/2(3 + p)(3_p)l/2 
128pl3/2(3+p)2(3_p)3/2 

X (-122553 - 15552p + 27270/ + 284V - 1593/ - 132/ + 4/) , 
while in the (2|0) sector we find 



pr^ = -A^, (A.11) 



(,o) , 1107 + 108. - 117/ -6/ 

1 8/(3+/ ' ^ ' 

p(2|0) ^ _>^3 ^ ^ 

^ 64p8(3+/^(3-p) 

X (-442341 - 65448p + 102330/ + 12924/ - 6669/ - 636/ + 80/) . (A.13) 

One of the main features of using multi-parameters transseries is the appearance of general- 
ized multi-instanton sectors, which may have different signs of the instanton action within the 
nonperturbative exponential contribution. This may lead, sometimes, to the cancellation of all 
terms in this exponential contribution — for example, in the present setting this happens when 
n = m — and we will be left with a (perturbative) expansion in the closed string coupling. The 
first sector with this feature is the (1|1) sector, where the first few coefficients are 

e) = A^, (A.14) 
_ x3 70713 - 10125p - 4617p^ + 261/ 

^2 - ^ ^^^8 ' {A.li>} 



J^l^) = (8806981203 - 1369011699. - 959100102/ + 103563198.^+ 

±.ZiO J) 

+18833715.^ - 787563/) . (A. 16) 



The general case n 7^ m is more complicated. Generically, asymptotic expansions will be in 
powers of the string coupling, Qs, and the "starting genus" may start taking negative values. 
Furthermore, logarithmic contributions begin to appear [15]. For example, one may compute 
the following coefficient in the (2|1) sector (this is the second non-vanishing coefficient in this 
sector): 

Pq^^'^^ = ^ ^ I (-432 - 180. + 24/ + 12/) + 

° 16//2(3+/(3-/'/H ^ ^ 

+ (153 - 15/ + 2/) log (^^^^ I . (A.17) 

Even though we have not produced as much data as in the one-cut solution discussed in 
[15], we are still able to conjecture the general form of all these coefficients. Our data, together 
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with the experience gathered in [15], indicate that the nonperturbative coefficients take the form 

min(n,m) 

pH-)(^)= ^ log'=(/(x))-P,("l™)W(x), (A.18) 

fc=0 

where 

p(n|m) [k] ( ) ^ >^ W (x) , (A. 19) 

and where the function f{x), written in terms of the variable p, is 

r,6 



Above, the coefficients Cj are given by 



ci = n + m + g - 1, (A.21) 
C2 = ^ (n + m) + 3^ - 1, (A.22) 

C3 = (1 - -Jnm) ^ (3m - n + 25) , (A.23) 
C4 = (1 - 6nm) (m + g), (A.24) 

and they are valid whenever n > m. The *Pg"'™'^^'^^ (x) are polynomials in p of degree 3 (m + g). 
When n = m, these polynomials get reduced and have degree n + g. Concerning the pattern for 
the coefficients, we find a similar result, but with the roles of C3 and C4 interchanged, 

\Cl 

and with the extra condition 

£j(n|m)[fc](^) = -qjC«l'^)W(_p). (A.26) 

Finally, upon further analyzing our data, a relation emerges between the coefficients in the 
(n|m)[A;] and the {n — k\m — k)[0] sectors (this is very similar to the relation (5.26) which we have 
found for the nonperturbative Painleve II coefficients in the main body of the text). We find 

p(n|m)[fc] _ 1 X{n-m) Y {n-k\m-k)[0] 

^ ~ k\\ 6 ) ■ 

For completeness, let us be fully specific on a few of the polynomials sp^'^l'^-'I^l^ which we 
have explicitly computed. These polynomials take the form c^ajp*, where p is their variable 
and c and their coefficients. We list these coefficients in the tables that follow. 

B. Perturbative Free Energy in the Quartic Matrix Model 

In the main text we have discussed how the Euler-MacLaurin formula (suitably adapted to the 
period-two case) provides for a recipe in order to extract the genus g perturbative free energies, 
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9 





1 


2 


c 


3 
4 


1 

82 


3 
512 


pO 


36 


670680 


811753164 


pi 


15 


8991 


163196127 


p2 


-2 


-159732 


-277986654 


p3 


-1 


486 


-53318331 






11340 


33149088 


p5 




-261 


6037173 


p6 




-208 


-1519542 


7 
pi 






-264789 


P' 






16680 


Q 






3068 



9 


-1 

— 1 





-1 

1 


n 


1 
6 


1 

16 


3 
256 


/ 


1 


153 


122553 


pi 







15552 


P^ 




-15 


-27270 


p3 




2 


-2844 


P' 






1593 


p5 






132 


p6 






-4 



Table 4: Prefactor c and coefficients of the polynomials ^g^'^^'*^' (left) and ^g^'^^'^' (right). 
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c 


1 

4 


1 
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3 
128 


p' 


297 


528525 


2176342749 


p' 


54 


69255 


368793810 


/ 


-27 


-118584 


-766103913 


/ 


-4 


-12744 


-118696752 


p' 




7533 


96370155 


P' 




549 


13149378 


p6 




-94 


-4972095 








-549324 


pS 






81696 








5288 



9 


-1 
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c 


1 
6 


1 

8 


3 

64 


pO 


1 


369 


442341 


p' 




36 


65448 


p' 




-39 


-102330 


p3 




-2 


-12924 


4 

p 






6669 


p5 






636 


P' 






-80 



Table 5: Prefactor c and coefficients of the polynomials *Pg'^'^'"'^' 



(left) and qjf 



(right). 
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3 
16 


3 

128 


9 

2048 


pO 


351 


2998377 


14430217473 


p' 


54 


451980 


2674565406 


p' 


-33 


-654156 


-4961854665 


p3 


-4 


-82620 


-841785048 


4 

p 




39393 


605366703 


p5 




3600 


90717246 


/ 




-382 


-29788263 


7 

p 






-3633804 


p8 






442368 


Q 
P^ 






30856 



9 


-1 
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1 

8 


3 

64 


9 

1024 


pO 


1 


1269 


3852765 


pi 




126 


593892 


p' 




-135 


-876582 


p3 




-8 


-115236 


4 

p 






55341 


p5 






5544 


p6 






-572 



Table 6: Prefactor c and coefficients of the polynomials ^Pg'''^^'^' 



(left) and ^f^''>^'^ 



(right). 



J-g^^^\ out of the recursion coefficients in the orthogonal polynomial framework. However, it 
is important to notice that this method is computationally very time consuming (even more 
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so that in the one-cut case addressed in [15]) and thus ends up providing for less data in the 
resurgence tests than directly using the coefficients pj"'^^ or Q^^^^\ Nonetheless, we explicitly 
need to know these coefficients as they are used to determine the Stokes coefficient out of the 
large-order sequence (4.42). In here, we shall explicitly list a few of these results for the J-g^^^^ 
(as usual, written in terms of the variable p). We find 



-^0 



(0|0) _ 



576A2 



logf^j, (B.l) 



J-f 1°) = , X (B.3) 

^ 320p6(9_p2)2 

X (787320 - 174960p - 2150551)2 + 43740^^ + 18630/ - 3780p^ - 47l/) , 



4'^'' = .... ... X (B.4) 

„2 , 1 innnc;o/i/iic^„3 



1792 (9 _p2 

X (1214950653504 - 234633327264p - 653277037896p^ + 119905844184p-'+ 



+141553030437p^ - 24374010024p^ - 15592951332/ + 2467933272p^ + 
+895852062p^ - 125778744/ - 23861628p^° + 2857680^^^ + 181989^^^) . 

As the genus increases, the expressions become exponentially longer and we shall not show any 
more explicit formulae. However, our results do indicate a clear pattern for the perturbative 
genus g free energies: for genus g >2 they have the form 

j-(o|o)(A,p) = (B.5) 

^ ^ ' p6(g-l) (9 _p2)2(3-l) "9^-^^' ^ ' 

where ^g{p) is a polynomial in p of degree 6{g — 1). Finally, as we have discussed in subsection 
5.3, applying the double-scaling limit to these results, and taking two derivatives, yields a precise 
match with (the square of) the perturbative data arising within the Painleve II equation. 

C. The Painleve II Equation: Structural Data 

In this appendix we present some explicit results concerning the two-parameters transseries 
solution to the Painleve II equation. Let us recall that in subsection 5.1 we have solved this 
equation, (5.8), by introducing the ansatz 



+00 +00 

u{w,aua2) = E E al^a2™e-("-™)^/-' <^h„)(^«), (C.l) 

71=0 m=0 



with 

min(n,m) +00 

^in\m)H^ log'^^^-^^"!"^)^^^^. (C.2) 

A:=0 g=0 
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For shortness we introduce 

+00 

g=0 

As we discussed in the main text, this ansatz turns the original differential equation into a 
recursive equation for the coefficients u^"''"^''^'. In table 7 we show the maximum order in w to 
which we have calculated these coefficients. We have only listed the n > m cases, but we shall 
see below how the (m|n) and {n\m) sectors are trivially related. We shall also see that there is a 
relation between the [k]th. and the [0]th logarithmic sectors. In the following we will reproduce 
a few examples concerning all this data. 



m\ 
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2 


3 


4 





1000 


1000 


1000 


500 


500 


1 




1000 


100 


100 


100 


2 






100 


100 


100 


3 








100 


100 


4 










100 



Table 7: Maximum order in w for which wc have calculated Wg"'™'"'^'. 

The first few (n|0) sectors we found are: 

^[0] , 1 4 73 8 10657 12 13912277 /r. .^ 

$ ^Ls =1 w w , (C.4) 

16 512 8192 524288 ' ^ ^ 

^roi 17 , 1513 5 850193 7 407117521 g 

^(110)=--^^ + 18432^ " 5308416^ + 2038431744^ ^^'"'^ 

<o] 1 2 41 4 5461 1734407 o 925779217 m x 

The lowest ^f'\.^ are 

(n|l) 

^[0] o 2 291 6 447441 886660431 14 ,^ 

^ = -3-0; w^^ u;^* , (C.7) 

128 32768 4194304 ' ^ ' 

<o] 3 115 5 30931 7 4879063 „ 

$ -\ -\ . (C.8 

48 18432 663552 ^ ' 

The first time we encounter logarithmic terms is for n = 2, m = 1, where we have 

17 4 1513 fi 850193 o ,^ , 

^ 12 2304 663552 ^ ^ 

From the full list of data we computed, one finds a relation between the coefficients in sectors 
(n|m)[A;] (logarithmic) and (n — k\m — A:)[0] (non-logarithmic), which is the following 

^{n\m)[k] = 1 (8 {m - n) )^-^^"-fc|™-fc)[0]. (c.io) 

Finally, the sectors with n < m are very closely related to the ones with n > m as^^ 

for n > m. (Cll) 



^{m\n)[k] 



(n|m)[fc] 
"9 



^^Similarly to what was found in [15] for the Painleve I equation, we may suspect that this is just an "apparent" 
relation, only to be falsified at some high n, m and g (in [15] one had to go to n = 3, m = 4 and genus g = 11 to 
falsify it). However, all the data we have produced is consistent with this relation. 
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As a final note, we add that all our off-criticality transseries results, i.e., the results for the 
two-cuts quartic matrix model partially presented in appendix A, match the present transseries 
solution of Painleve II, when in the double-scaling limit. A Mathematica notebook with the 
complete explicit results we have obtained is available upon request. 
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